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The  thermal  conductivity  and  specific  heat  of/  NdP^O^  have  been  measured 
using  a dynamic  electrical  heating  technique.  J The  values  of  the  thermal  con- 
ductivity tensorelements  (in  W/cm-K)  are /aS-  4il0  J(  16”*,  = 9.66  x id^ 

and  iA^  = 1.40  x 10*2.  The  specific  heat  is  0Tl41  cal/g-K,  in  excellent  agree- 
ment with  the  value  of  0.139  cal/g-K  that  was  measured  by  differential  scanning 
calorimetry. J C = . <2>  I*/ 

OstrtgThe experimentally  mennneed -thermal  properties.  we  have  calculated  the 
time  -dependent  heating  effects  of  typical  flash-lamp-excited  round  and  rectan- 
gular laser  rods.  With  repetition  rates  limited  to  flO  ppm,  simple  conductive 
cooling  by  an  ambient-temperature  heat-sink  should  be  adequate. 

Computer  simulation  of  the  flash-lamp-excited  NPP  laser  has  been  carried  out 
by  numerical  solution  of  the  laser  rate  equations.  With  1 J electrical  input, 
typical  calculated  output  energies  are  24  mj  for  normal-mode  operation  and 
10  m J for  passively  Q-switched  operation.  Optimum  Q-switch  absorption  and 
output  mirror  transmission  have  also  been  calculated. 

Experimental  studies  of  such  lasers  have  been  made  with  varying  laser  rod 
sixe  and  laser  resonator  parameters.^- With  the  largest  rod,  2 x 2 x 20  mm\ 
we  have  obtained  72  mj  output  with  jTw  input,  in  a normal-mode  confocal 
resonator.  Passive  Q-switchlng  of  2 x 2 x 10-mm^  rods  in  a hybrid  unstable 
resonator  has  produced  5 mj  near-field  output  (2  J input)  in  a 6-nsec  pulse, 
with  3.9  mJ  collected  in  a 3.1-mrad  far-field  aperture. 
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I 

I.  INTRODUCTION 

This  report  describes  the  results  of  research  studies  on  neodymium  pentaphosphate 
minature  lasers.  The  high  Nd  concentration  in  this  material  (NdP,.014,  abbreviated  as  NPP) 
allows  effecient  absorption  of  pump  radiation  and  the  construction  of  small  laser  devices  which 
should  be  suitable  for  rangefinder  applications.  Most  of  the  important  fluorescence  and  optical 
properties  of  NPP  have  been  previously  measured  and  reported,  so  that  our  efforts  were  con- 
fined to  two  other  general  areas:  the  determination  of  the  material's  thermal  properties  and 
the  evaluation  of  flash-lamp-excited  laser  performance.  In  Section  II  of  this  report,  we  de- 
scribe the  method  and  results  of  the  measurement  of  the  thermal  conductivity  and  specific  heat 
of  NPP.  Using  these  data,  we  also  do  some  simple  calculations  to  estimate  the  thermal  limi- 
tations of  NPP  laser  performance.  In  Section  III,  we  examine  the  characteristics  of  NPP  flash- 
lamp-excited  lasers.  First,  the  details  of  a computer  simulation  of  such  lasers  are  given,  in- 
eluding  the  effects  of  varying  the  Nd  concentration  in  the  laser  rod,  and  of  incorporating  a 
saturable -absorber  Q-switch  in  the  laser  cavity.  Next,  the  experimental  results  describing 
the  operation  of  NPP  lasers  are  presented.  These  results  include  operation  in  normal  and 
Q-switched  modes,  ip  stable  and  unstable  resonators,  and  with  varying  laser  rod  sizes.  . Finally, 
in  Section  IV  we  summarize  the  results,  and  indicate  further  areas  of  investigation  of  assess- 
ment of  NPP  as  a rangefinder  laser  material. 

All  of  the  NdPgOj^  crystals  and  laser  rods  described  in  this  report  were  provided  by 
Dr.  W.  K.  Zwicker,  Philips  Laboratories,  Briarcliff  Manor,  New  York. 


II.  THERMAL  PROPERTIES 


One  of  the  problems  in  operation  of  high-energy,  fast-repetition-rate  pulsed  lasers  is  heat 
removal  from  the  laser  rod.  In  order  to  analyze  this  problem,  we  have  measured  the  thermal 
conductivity  and  specific  heat  of  NPP  using  a recently  published  technique.  Most  of  the  perti- 
nent analysis  and  results  are  given  in  the  following  Section  II-A,  parts  (1)  - (5)  of  which  are  the 
manuscript  of  an  article  (with  W.  K.  Zwicker)  to  be  published  in  the  Journal  of  Applied  Physics. 
Additional  experimental  details  and  data  are  found  after  the  first  part,  and  an  elaboration  of  the 
theoretical  calculations  of  the  conductivity  measurement  is  found  in  the  Appendix.  In  Section  R-B, 
we  present  some  calculations  on  laser  rod  heating  based  on  the  measured  thermal  conductivity. 

A.  THERMAL  CONDUCTIVITY  AND  SPECIFIC  HEAT  OF  NdP,0,,. 

5 14 

1.  Background 

Neodymium  pentaphosphate  is  a high -Nd -concentration  material  with  small  enough  fluores- 
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cence  quenching  for  use  in  miniature,  efficient  lasers.  ' Growth  of  crystals  of  NPP  larger  than 
previously  available  and  having  volumes  of  several  cubic  centimeters3  has  permitted  the  fabrica- 
tion of  rods  for  flash -lamp-excited  lasers.4  One  objective  of  the  further  development  of  high-Nd- 
concentration  lasers  is  to  make  CW  devices,  pumped  either  by  lamps  or  by  light -emitting  diodes. 
Another  goal  is  to  increase  the  repetition  rate  and  average  power  of  the  flash-lamp-pumped 
lasers.  For  either  of  these  developments,  heat  removal  from  the  laser  rod  is  of  prime  concern, 
and  the  thermal  conductivity  of  the  rod  must  be  known  for  proper  design.  We  have  measured  the 
anisotropic  thermal  conductivity  tensor  and  the  specific  heat  of  NPP  by  using  a simple  extension 
of  the  technique  recently  described  by  Bruce  and  Cannell5  for  determining  the  thermal  conduc- 
tivity and  specific  heat  of  isotropic,  electrically  insulating  solids.  The  basic  principle  of  mea- 
surement is  straightforward.  Two  parallel  metal  strips  are  fabricated  on  the  surface  of  the 
sample.  Alternating  current  is  passed  through  one  strip,  the  heater,  and  the  temperature- 
induced  periodic  resistance  variation  of  the  other  strip,  the  sensor,  is  measured  as  a function 
of  the  heating  frequency.  From  the  measured  data,  and  together  with  the  electrical  parameters 
for  the  heater  and  sensor  and  the  temperature  coefficient  of  sensor  resistance,  the  thermal  con- 
ductivity and  specific  heat  can  be  found. 

The  theoretical  relation  of  the  thermal  conductivity  tensor  elements  to  the  measured  data 
is  given  in  Section  II-A-2  below.  In  Section  II-A-3,  we  briefly  describe  the  experimental  method, 
particularly  the  aspects  in  which  it  differs  from  the  method  of  Ref.  5.  The  experimental  results 
are  presented  in  Section  II-A-4  and  discussed  in  Section  II-A-5. 

2.  Theory 

We  have  extended  the  theoretical  calculation  of  Ref.  5 to  allow  for  different  thermal  conduc- 
tivities along  the  different  crystallographic  axes.  A related  derivation  for  transient  excitation 
is  found  in  Ref.  6.  Since  the  space  group  of  NPP  is  monoclinic,7  in  principle  it  is  necessary  to 
find  a new  set  of  principal  axes  along  which  the  thermal  conductivity  tensor  is  diagonal.8  How- 
ever, the  monoclinic  angle  0 deviates  from  90*  by  only  0.S*.  so  we  shall  make  the  simplifica- 
tion of  assuming  orthorhombic  symmetry,  with  the  tensor  diagonal  along  the  three  nearly  per- 
pendicular crystal  axes,  and  having  components  A{  = Aft,  Ab>  and  Ac- 
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The  following  discussion  uses  the  notation  of  Ref.  5 as  closely  as  possible,  and  proceeds 
along  similar  lines,  except  that  we  first  find  the  slightly  more  general  three-dimensional,  time 
dependent  Green's  function  (see,  for  example.  Ref.  9).  The  heat-flow  equation  becomes 


T = temperature 
p = density 

C = specific  heat  (heat  capacity/unit  weight) 
P(x,t)  = power  dissipation/unit  volume. 

The  equation  for  the  Green's  function  is 


4»4(r  — r ) «(t-t  ) 


with  solution 


[The  function  u(t)  is  the  caus  .l  unit  step  function.]  With  a power  excitation  density  of  PQ4(r) 
cos  ait  in  Eq.  (1),  integration  over  time  gives  the  solution  (see  Appendix) 

P „/  2 2 2 \-l/2  . „ 

+ b + I-)  (KKK)'i/Z 


(Note  that  u,  the  power  radian  frequency,  is  at  twice  the  frequency  of  the  actual  heater  current.) 
We  now  assume,  as  in  Ref.  5,  that  the  heat  source  distribution  (coming  from  the  heater  strip)  is 
a line  in  the  plane  at  z = 0 which  extends  in  the  y direction  from  y = — « to  +*°.  The  temperature 
in  the  z = 0 plane,  at  a distance  d in  the  x direction  from  the  line  source  (corresponding  to  the 
sensor  strip  location)  is  found  by  integrating  Eq.  (4)  over  the  y coordinate.  The  temperature 
at  the  sensor  is 


where  Kq  is  a modified  Bessel  function  of  imaginary  argument,  and  P is  the  heater  power/unit 
length.  In  Eq.  (5),  we  have  also  taken  the  solution  for  the  case  of  a heater  source  on  the  surface 
of  the  infinite  half-volume  z < 0,  using  the  same  symmetry  arguments  as  in  Ref.  5. 

The  Important  point  to  be  noted  is  that  Eq.  (5)  is  the  same  as  Eq.  (6)  of  Ref.  5 except  that 
<AxAzf 1/2  has  replaced  the  initial  factor  A-1  and  (aipC/A^1/2  has  replaced  (wpC/A)1//2  in  the 
argument  of  the  Bessel  function.  These  substitutions  are  physically  reasonable  because  the 
dynamic  response  of  the  heater-sensor  pair  given  by  the  frequency-dependent  Bessel  function 
depends  only  on  the  dlffuslvity  in  the  direction  between  them,  whereas  the  magnitude  of  the 


response  given  by  the  initial  factor  depends  on  the  relative  amounts  of  heat  flow  in  both  x and 
z directions  out  from  the  heater  strip.  Thus,  the  smaller  the  downward  conductivity,  Az>  the 
larger  will  be  the  heat  flow  and  response  in  the  x direction.  Of  course,  the  A^  component  does 
not  appear  since  with  the  assumed  symmetry  and  infinite  y extent,  there  is  no  heat  flow  in  the 
y direction. 

These  arguments  also  apply  to  the  correction  factor  discussed  in  Ref.  5 for  the  non-zero 
width  of  the  heater  and  sensor  strips,  so  that  the  final  frequency-dependent  complex  voltage 
which  is  measured  across  the  sensor  is 


where 


V = 


VdP^  KjQd)  Isinh|ow/2)]2 


V = DC  bias  voltage  across  sensor 

/?  = temperature  coefficient  of  sensor  resistance,  (1/R)  dR/dT 
a = (iu)pC/Ax)1'/Z 

d = separation  between  heater  and  sensor  strips 
w = width  of  strips. 


The  magnitude  of  the  voltage  response  is  similar  to  Eq.  (15)  of  Ref.  5,  that  is: 


VPfl  j^sinh2  (aw/2  \TZ)  + sin2  (aw/2  v/T)|  |K  (Q<j|)| 
<5w/2>2  ' °“ 


(6) 


(7) 


where  a = (upC/A^)  ' and  the  quantity  A has  become  J 

For  various  crystal  orientations  with  respect  to  the  heater  and  sensor  strips,  the  appro- 
priate components  A^  = A&,  A^,  A are  used  in  the  above  equations,  assuming  the  simple  geom- 
etry of  principal  axes  oriented  parallel  and  perpendicular  to  the  strips. 

3.  Experimental  Procedure 

The  measurement  apparatus  differed  only  in  minor  details  from  that  described  in  Ref.  5. 
For  example,  we  used  an  additional  two-phase  lock-in  amplifier  to  monitor  the  heater  current. 

The  sensor  temperature  coefficient  p was  found  by  heating  each  sample  and  measuring  R 
over  the  range  from  23°  to  50° C.  Polynomial  regression  was  used  to  obtain  the  best  fit  to  the 
data  and  to  determine  p = (l/R)  dR/dT  at  the  temperature  of  interest,  27  °C. 

The  material  used  for  heater  and  sensor  strips  was  sputtered  Ni,  500  A thick.  The  strip 
patterns  were  defined  by  conventional  masking  and  photolithographic  techniques,  and  the  un- 
wanted Ni  areas  were  removed  either  by  chemical  etching  or  film-parting.  A layer  of  Au  was 
evaporated  over  the  contact  pads  at  the  ends  of  the  two  strips,  and  0.001 -in. -dia.  Au  bonding 
wires  were  attached  to  the  contact  pads,  using  acoustic  bonding  with  silver-epoxy  support,  or 
just  silver-epoxy  alone.  In  all  cases,  the  ratio  of  heater-to-sensor  resistance  was  the  same 
as  the  ratio  of  heater-to-sensor  length,  showing  that  the  external  connection  resistances  were 
negligible.  Figure  1 is  a photomicrograph  of  a typical  str*p  pattern,  with  the  external  bonding 
wires  removed.  The  strip  patterns  were  generally  smaller  than  those  of  Ref.  5,  with  minimum 
lengths  of  2.5  mm,  widths  of  ~4  pm,  and  separations  ~33  pm. 
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Fig.  1.  Photomicrograph  of  heater  and  sensor  strips  with  bonding  pads, 
on  a cleaved  b face  of  NPP  (the  bonding  wires  have  been  removed).  The 
heater  (top)  is  3 mm  long,  and  the  sensor  (bottom)is  2.5mm  long.  The 
strips  are  4 jim  wide,  33  |im  apart,  and  run  parallel  to  the  a direction. 


All  samples  were  cut  or  cleaved  from  the  same  NPP  crystal.  To  determine  the  three  ther- 
mal conductivity  tensor  elements  and  the  specific  heat,  three  different  orientations  were  used. 
The  results  for  samples  of  the  same  orientation  agreed  within  a few  percent,  even  with  consider 
able  variation  in  strip  parameters  and  power  dissipation.  The  thermal  properties  were  derived 
from  the  results  for  the  one  sample  of  each  orientation  that  gave  voltage  versus  frequency  data 
showing  the  smallest  deviation  from  the  best-fit  theoretical  curves.  Sample  Ml  had  a polished 
c face  with  the  metal  strips  in  the  a(y)  direction.  In  this  case,  Ax  = A^  and  Az  = Ac>  Sample 
M2  also  had  a polished  c face,  but  with  orthogonal  strip  orientation,  i.e.,  Ax  = A&  and  Az  = 
Sample  M3  had  a cleaved,  step-free  b face  with  strips  in  the  a direction,  giving  A = A and 


4.  Results 


For  each  sample,  data  at  27*C  were  obtained  over  a range  of  roughly  two  orders  of  magni- 
tude in  both  the  heater  frequency  and  sensor  voltage  response.  A best  fit  to  the  data  was  ob- 
tained by  varying  the  two  independent  parameters  and  pc/ Ax  to  minimize  the  least- 

squares  difference  between  the  logarithm  of  the  measured  values  of  | | and  the  logarithm  of 

the  right  side  of  the  theoretical  expression,  Eq.  (7).  (This  was  done  to  ensure  comparable, 
although  smaller,  weighting  of  the  high-frequency  points,  where  the  absolute  voltage  response 
is  smallest.)  No  weighting  was  carried  out  to  compensate  for  frequency-dependent  noise. 4  5 Typ- 
ically, the  largest  deviation  between  the  measured  and  calculated  values  of  | V | was  ~1.5  per- 
cent. A set  of  representative  data  (for  M3)  is  shown  in  Fig.  2,  with  its  best-fit  curve.  The  peak 
value  of  | Vw  | is  plotted  as  a function  of  f = u>/2ir,  the  heating  frequency  (twice  the  current  drive 
frequency).  For  this  particular  set  of  data,  the  experimental  parameters  were 

V = 0.869  V 
P = 2.07  mW/cm 


f -w/2*  (M«» 


Fig.  2.  .Sqnpor  voltage  response  | V | as  a function  of  heater  frequency 
f = o)/2ir  for  NPP  sample  M3.  Experimental  parameters  are  given  in 
the  text.  The  open  circles  are  the  data  points,  and  the  solid  curve  is 
calculated  using  the  indicated  best-fit  parameters  for  ^ AfeA  and  pC/A  . 


The  two  dependent  parameters,  J~K~K~  and  pC/A^,  obtained  for  each  of  the  three  orientations, 
are: 


Sample 


./A^yW/cm-K) 


Ml  = 1172  * 10 


M2  = 1-713  x 10 

M3  = 1-159  x 10 


(dlffusivity)-1  (sec-cm~2) 


pC/Ab  = 204.9 
pC/Aa  = 94.3 
pC/Ac  = 141.2 


The  uncertainty  in  the  values  for  NyrA~A’  is  about  ±2  percent,  and  arises  from  the  cumulative 
errors  in  measuring  | Vw  | , P,  and  0.  The  diffusivities,  A^/pC,  are  independent  of  the  absolute 
system  response  calibration  factors.  We  estimate  that  the  diffusivities  are  accurate  to  about 
±1  percent,  reflecting  the  goodness  of  fit  to  the  shape  of  the  | V (f)|  curve.  (On  the  logarithmic 
plot,  the  initial  multiplicative  factor  shifts  the  curve  without  changing  its  shape.  ) The  thermal 
conductivity  ratios  A^/  Aj  can  be  found  from  the  diffusivities  without  knowledge  of  C or  absolute 
system  response  factor. 

The  above  data  are  redundant  in  that  J~A^A~c  is  obtained  for  two  orientations;  the  two  values 
agree  to  within  1 percent,  confirming  the  overall  calibration  accuracy.  For  the  further  data  re- 
duction, we  use  the  average  value  J \AC  = 1.166  X 10"2  W/cm-K.  The  value  of  Aa/Ab  = 2.173 
found  by  taking  the  ratio  of  diffusivities  agrees  to  better  than  1 percent  with  value  Afi/Ab  = 2.15P 
found  by  taking  the  square  of  the  ratio  of  to  J A^  A c from  the  A factors. 

Combining  the  J A^Ac  and  diffusivity  results  for  samples  Ml  and  M3,  we  can  calculate  the 
specific  heat  (per  unit  volume),  pc,  as  follows: 


AbAc  = (1.66  X 10‘2)2  = (pC/ 204.9)  x (pC/141.2)  . 

giving  pC  = 1.98  j/cm3-K.  Similarly,  for  M2  and  M3. 

AaAc  = (1.713  X 10"2)2  = (pC/94.3  ) X (pC/141.2)  , 

also  giving  pC  = 1.98  j/cm3-K.  The  high -temperature  Dulong-Petit  limiting  value  for  pC  in 
NPP  is  3.21  j/cm3-K.  Using  the  measured  value  of  pC  with  the  diffusivity  parameters  gives 
the  following  values  of  the  thermal  conductivity  tensor  elements: 

Aa  = 2.10  x 10’2  Afa  = 9.66  X 10‘3  Ac  = 1.40  X 10-2  W/cm-K  . 

The  density  of  NPP  calculated  from  its  lattice  parameters^  is  3.365  g/cm3,  which  gives  the 
specific  heat  (per  unit  weight)  (£  = 0.141  cal/g-K. 

Using  small  pieces  from  the  same  NPP  single  crystal  as  above,  we  also  measured  the 
specific  heat  directly  with  a DuPont  900  Differential  Thermal  Analyzer  equipped  with  a differ- 
ential scanning  calorimetry  (DSC)  cell.  Powders  of  Zn  and  A1  were  used  as  specific-heat  cali- 
bration materials.  The  directly  measured  specific  heat  of  NPP  was  <?  = 0.139  ± 0.002  cal/g-K, 
in  excellent  agreement  with  the  value  obtained  by  the  electrical  dynamic  measurement. 

5.  Discussion 

The  thermal  conductivity  of  NPP  is  low  compared  with  that  of  YAG,  one  of  the  better  heat- 
conducting  laser  hosts.  Even  the  highest  component  in  NPP,  Aft  = 2.10  x 10'2  W/cm-K,  is  only 
one-sixth  of  the  value  for  YAG,  0.13  W/cm-K  (Ref.  10).  However,  this  does  not  rule  out  NPP 
for  CW  or  high-repetltlon-rate  lasers,  since  rods  of  NPP  will  typically  be  much  smaller  than 
their  Nd:YAG  counterparts.  Even  in  very  small  Nd:YAG  laser  systems,  the  rod  diameter  is 
commonly  at  least  3 mm  to  ensure  efficient  absorption  of  pump  radiation.  Since  NPP  has 
30  times  higher  Nd  concentration  than  NdrYAG  with  standard  1%  Nd  doping,  smaller  NPP  rod 
diameters  can  be  used.  For  example,  a 1-mm2  rod  was  used  in  Ref.  4.  The  temperature  rise 
at  the  rod  center,  for  a given  heat-sink  surface  temperature,  varies  as  the  square  of  the  rod 
dimension  and  inversely  with  the  thermal  conductivity.11  By  using  a NPP  rod  with  about 
1 /V?  ~ 0.4  times  the  transverse  dimension  of  a Nd:YAG  rod,  comparable  cooling  can  be 
achieved. 

The  thermal  conductivity  of  NPP  is  rather  anisotropic,  with  A :A  :A^  = 2.17:1.45:1.  This 
thermal  anisotropy  is  consistent  with  both  the  mechanical  properties  of  NPP  crystals  and  the 
microscopic  crystal  structure.  The  lowest  conductivity  direction  is  perpendicular  to  the  easy- 
cleavage  b plane,  and  the  highest  conductivity  is  in  the  a,  direction,  along  the  structurally 
close-knit  ribbons  of  phosphate  groups.  If  cooling  of  the  laser  rod  is  nonuniform,  over  only 

part  of  the  rod  surface,  it  is  therefore  important  to  ensure  that  the  heat  flow  is  in  the  a 
direction. 

The  present  technique  for  electrical  dynamic  measurement  of  the  thermal  conductivity  and 
specific  heat  Is  quite  simple  and  convenient.  For  measurement  of  these  properties  in  NPP,  in 
particular,  it  afforded  several  advantages: 

(a)  The  sample  size  could  be  quite  small,  with  measurement  surfaces  about 
2 mm  X 4 mm  and  a minimum  thickness  of  only  about  2 mm. 


8 


r 


su:; 


(b)  All  three  thermal  conductivity  tensor  elements  could  be  determined 
using  only  the  two  b and  c crystal  surfaces  for  measurement.  In  a 
conventional  direct  heat  flow  measurement,  it  is  desirable  to  have 
smooth,  flat  sample  surfaces  to  minimize  thermal  impedance  across 
the  interfaces.  Measurement  of  Aa  by  the  direct  technique  would  have 
required  a polished  a surface,  which  is  difficult  to  prepare. 

(c)  The  specific  heat  was  determined  along  with  A.,  providing  a parameter 
which  could  be  compared  with  the  calorimetrically  measured  value  for 
confirmation  of  the  experimental  accuracy. 

The  transient  method  of  Ref.  6 would  also  offer  these  advantages,  but  experimentally  it  would 
probably  be  more  difficult  to  obtain  accurate  measurements  of  the  low-level  time -dependent 
voltage  response  than  the  steady-state  frequency-dependent  response.  An  interesting  technique 
for  determining  the  thermal  diffusivity  by  measuring  laser-induced  temperature  distributions12 
might  also  be  applicable,  but  would  be  complicated  by  the  fluorescence  properties  of  NPP. 

6.  Additional  Experimental  Details  and  Data 

Block  diagrams  of  the  experimental  apparatus  are  shown  in  Figs.  3 and  4.  In  Fig.  3 is  given 
the  setup  for  determining  the  sensor  temperature  coefficient.  The  temperature  controller  sta- 
bilizes the  sample  (and  sensor)  temperature,  which  is  measured  by  a calibrated  Pt  resistance 
and  digital  ohmmeter.  The  resistance  of  the  sensor  strip  is  measured  using  a DC  current 
source  set  at  100  ^A  and  measuring  the  voltage  with  the  digital  voltmeter.  The  resistance 
values  were  measured  from  room  temperature  to  ~50*C. 

The  apparatus  for  measuring  the  thermal  conductivity  is  shown  in  Fig.  4.  The  low-frequency 
function  generator  provides  the  AC  power  to  the  heater  resistance  element,  Ru.  A calibrated 
iOn  resistance  in  series  with  the  heater  is  used  to  determine  the  current  through  R^,  and  hence 
the  heating  power.  A battery-powered  current  source  provides  the  DC  bias  to  the  sensor  ele- 
ment, Rg,  and  the  AC  temperature-induced  voltage  component  is  amplified  and  detected  at  twice 
the  heater  frequency  by  the  lock-in  amplifier. 

The  experimental  parameters  for  the  samples  Ml  and  M2,  not  given  above,  are: 


Ml 

M2 

V 

0.681 

1.404 

V 

p 

1.73 

3.17 

mW/c: 

0 

1.80  x 10‘3 

1.35  X 10"3 

K'1 

w 

8.19  X 10'4 

4.68  X 10'4 

cm 

d 

6.35  X 10"3 

3.27  X 10'3 

cm 

The  experimental  data  for  Ml  and  M2  are  shown  in  Figs.  5 and  6,  with  the  best-fit  curves 
and  parameters  and  pC/A.  for  each  case. 
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Fig.  3.  Experimental  apparatus  for  measuring  temperature  coefficient 
of  sensor  element  on  NPP  sample. 


Fig.  4.  Experimental  apparatus  for  measuring  thermal  conductivity 
and  specific  heat  of  NPP. 
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Fig.  5.  Sensor  voltage  response  as  a 
function  of  heater  frequency  for  NPP 
sample  Ml. 


Fig.  6.  Sensor  voltage  response  as  a 
function  of  heater  frequency  for  NPP 
sample  M2. 


B.  TIME-DEPENDENT  HEATING  EFFECTS 

1.  Statement  of  the  Problem 

The  problem  that  is  addressed  in  this  section  is  the  determination  of  the  temperature  rise 
in  a pulsed  NPP  laser,  and  the  resulting  limitations  on  pump  energy  or  repetition  rate.  In 
principle,  such  calculations  can  be  done  as  accurately  as  desired,  but  require  a good  knowledge 
of  the  excitation  geometry  and  rod  cooling  parameters,  and  a subsequent  numerical  solution  of 
the  heat-diffusion  equation.  For  the  present,  we  are  primarily  concerned  with  an  easily  ob- 
tained semi -quantitative  answer  to  the  heating  problem,  so  that  we  restrict  ourselves  to  simple 
models  for  which  analytic  solutions  can  be  found. 

First,  the  experimental  situations  which  we  have  dealt  with,  and  presumably  which  would 
be  encountered  in  the  laser  rangefinder,  are  for  uniform  transverse  mode  structure  across  the 
laser  rod.  This  could  arise  either  from  multimode  alignment  of  a stable  resonator,  or  a simple 
uniform  transverse  mode  of  an  unstable  resonator.  For  this  reason,  we  choose  the  initial  ther- 
mal excitation  to  be  uniformly  distributed  over  the  rod  cross  section.  This  is  probably  not  the 
case,  due  to  nonur.iform  absorption  of  the  pump  radiation,  but  even  a nonuniform  temperature 
distribution  would  have  to  be  integrated  over  the  laser  mode  profile  to  give  the  total  average 
heating  effect.  Also,  the  details  of  finding  the  nonuniform  excitation  profile  and  of  then  calcu- 
lating the  time  evolution  would  be  quite  complicated. 

Second,  we  approximate  the  flash-lamp  excitation  to  have  a delta-function  time  dependence. 
This  is  probably  a good  approximation  because,  as  we  shall  see,  the  lamp  pulse  duration, 
~50  (jisec,  is  very  much  shorter  than  any  characteristic  heat  diffusion  time  in  these  lasers. 

Third,  we  shall  treat  only  simple  rod  geometries.  First,  we  examine  the  cylindrical  rod 
with  isotropic  conductivity  and  uniform  heat-sinking.  Next,  we  treat  the  problem  of  one- 
dimensional heat  flow  across  a slab.  The  actual  case  of  a rectangular  rod,  heat  sunk  on  one 
surface  and  pumped  primarily  on  the  opposite  face,  represents  a combination  of  these  two  cases. 


11 


2.  Cylindrical  Rod 

Taking  the  Laplace  transform  of  Eq.  (1),  in  cylindrical  coordinates,  with  radial  symmetry, 
we  obtain 


where 


q = p(pC/A) 

For  a delta-function  excitation, 


with  the  form  of  solution' 


Since  v is  finite  at  r = 0,  c^  = 0.  The  other  boundary  condition  is  that  the  temperature  be  con 
stant  (here  set  to  zero)  at  the  boundary  r = a.  The  boundary  conditions  and  Eq.  (9)  give 


The  inverse  transform  is 


with  p = *J\pC/ A,  which  is  evaluated  by  finding  the  residues  of  the  contour  integral  in  the  left 
half  complex  plane.  The  result  is 


where  a are  the  zeros  of  the  Bessel  function  given  by  J (o  a)  = 0.  In  dimensionless  units 

n n n 


The  temperature  is  plotted  as  a function  of  normalized  spatial  coordinate  u In  Fig.  7,  at  fixed 
Intervals  from  t = 0.05  to  0.50  in  increments  of  0.05,  and  has  itself  been  normalized  by  the 


I 


Fig.  7.  Temperature  as  a function  of  radial 
distance  for  different  times  after  a uniformly 
distributed  heat  impulse  in  a cylindrical  rod. 
All  parameters  are  normalized  as  explained 
in  the  text. 


Fig.  8.  Logarithmic  dependence  of  tempera- 
ture at  rod  center  with  time  (solid  circles). 
The  time-dependence  for  the  rectangular  slab 
is  also  shown  (open  circles),  with  both  limit- 
ing decay  times  being  indicated. 


parameter  (b/pC).  After  a short  interval  (T  > 0.2),  the  temperature  at  the  rod  center  is  domi- 
nated by  the  lowest  order  normalized  decay  constant,  rQ  = l/u^  = 0.173,  so  that  this  single  ex- 
ponential decay  is  a good  approximation  to  the  subsequent  temperature  evolution,  as  shown  in 
Fig.  8. 

Taking  some  typical  parameters,  a = 1 mm  and  an  average  diffusivity  A/pC  = 7.5  x 
10  3 cm^/sec,  the  actual  decay  time  tq  = 0.23  sec.  The  initial  temperature  rise  is  b/pC  which 
can  be  estimated  from  the  pumping  parameters.  Of  the  1 J of  electrical  energy  into  the  flash- 
lamp,  approximately  60  percent  is  radiated,  and  we  will  assume  that  half  of  this  is  coupled 
into  the  laser  rod.  On  the  order  of  20  percent  of  this  broadband  radiation  is  absorbed  by  the 
Nd3+  ions,  giving  b »60  mj.  Approximately  10  mj  might  be  emitted  as  laser  radiation,  giving 
a net  value  of  B = 50  mJ  of  thermal  energy  input.  Since  pC  » 2 J/cm3-K,  the  temperature  rise 
is 


AT  = 


B 1 _ 0.05  J 1 A 0 „ 

T—  T7T  = J \ 1 = 0.8  K 

(»a‘i)  T(0.1  r X 1 cni  2 j/cm  -K 


In  a periodic  pulse  train  with  period  T.  with  simple  exponential  decay,  the  temperature  will  in- 
crease by  a factor  of  [1  - exp(-T/ro)I  \ which  should  be  negligible  for  T > 4tq  =»  1 sec.  Of 
course,  this  calculation  neglects  convective  heating  by  the  lamp,  or  other  types  of  radiation 
absorption  processes  in  the  rod  such  as  direct  UV,  atomic,  or  infrared  lattice  absorptions. 


3.  Rectangular  Rod 

In  examining  the  case  of  the  rectangular  laser  rod,  we  shall  simplify  the  situation  and  as- 
sume that  no  heat  flow  occurs  across  surfaces  which  are  not  contacted  to  the  heat  sink.  This 


13 


transforms  the  problem  to  a simple  one -dimensional  model,  and  we  will  examine  the  situations 
of  having  either  one  face,  or  two  opposite  faces  cooled. 

The  simple  rectangular-coordinate  heat  flow  equation,  with  no  y variation,  corresponding 
to  Eq.  (9 ) is 


with  solution  of  the  form 


v = c,  sinhqx  + c?  coshqx  + c 


First,  apply  the  boundary  condition  that  at  two  faces  the  temperature  is  constrained  to  be  zero 
v(x  = ±L)  = 0,  which  gives  the  solution 


again  with  p = VXpC/A.  Evaluation  of  the  contour  integral  by  summation  of  the  residues  results 


where 


At  this  point,  we  can  also  find  the  solution  for  the  slab  with  only  one  face  fastened  to  the  heat 
sink,  and  no  heat  flow  across  the  opposite  face.  It  is  convenient  to  shift  the  coordinate  system 
so  that  the  boundary  conditions  are  v(2L)  = 0 and  (8v/8x)  (0)  = 0.  With  these  constraints 


c.  sinh  2qL  + c,  cosh2qL  + c,  = 0 


c^q  cosh(0)  + c2q  sinh<0)  = 0 

giving 

V = pCp  fl  " ^shZqll  * (2°) 

As  might  be  expected  from  symmetry  arguments,  Eq.  (20)  and  its  inverse  transform  are  identical 
to  Eqs.  (17)  and  (19)  with  the  replacement  of  the  dimension  L by  2L,  or 


u = x/2L 


un  = (n  + 1/2)  * 


■y  A t 


The  time  evolution  of  the  normalized  temperature  is  shown  in  Fig.  9.  In  this  case,  the  asymp- 
totic normalized  decay  constant  is  to  = 1/u2  = (»/2)~2  (Fig.  8),  and  tq  = (pC/A)  (2L)2  (2/ir)2.  For 
a rod  2 mm  thick  (2L  = 0.2  cm),  with  pC/ Aft  = 94.3  sec/cm2,  the  decay  time  is  tq  = 1.53  sec. 


Fig.  9.  Temperature  as  a function  of  distance  (either  from  slab  center, 
or  uncooled  slab  face)  for  different  times  after  a uniformly  distributed 
heat  pulse  in  the  slab,  for  rectangular  rod  with  no  transverse  heat  flow. 
The  boundary  conditions  determine  the  normalization  parameters. 


The  initial  temperature  rise  is 


(2L)  l 


0.05  J 1 

2 1 1 — 

(0.2)  x 1 cm  2 J/cm  -K 


= 0.63  K 


Although  at  first  glance,  this  one -dimensional  problem  may  seem  like  an  extreme  simplifi- 
cation, it  probably  is  more  accurate  a model  than  the  uniformly  cooled  cylindrical  rod.  Irre- 
spective of  the  actual  energy  input  to  the  rod,  the  thermal  time  constant  will  be  on  the  order  of 
1 or  2 sec.  With  very  close  coupling  to  the  flash  lamp,  convective  heating  may  also  be  impor- 
tant, as  well  as  overall  conductive  heating  of  the  entire  laser  head  assembly.  Therefore,  unless 
a more  elaborate  cooling  method  is  used,  simple  conductive  cooling  of  one  rod  face  will  probably 
limit  the  rangefinder  repetition  rate  to  about  10  ppm. 
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III.  NPP  FLASH-LAMP-EXCITED  LASER  STUDIES 

A.  COMPUTER  SIMULATION 

1.  Description  of  the  Problem 

i 

In  an  effort  to  improve  our  understanding  of  the  operation  of  flash-lamp-excited  NPP  lasers, 
we  have  carried  out  numerical  analysis  of  these  devices  by  means  of  computer  solution  of  the 
laser  rate  equations.  Such  simulation  should  be  useful  for  several  reasons: 

(a)  We  should  obtain  a semi -quantitative  estimate  of  expected  device 
efficiencies  under  different  operating  conditions  to  compare  with  our 
experimental  results. 

(b)  We  should  have  a guide  for  experimental  studies  and  for  optimizing  the 
actual  laser  parameters  by  studying  the  simulated  effects  of  varying 
Q-switch  saturable  absorber  and  output-coupling  parameters. 

(c)  We  can  calculate  the  effects  of  reducing  the  Nd  concentration,  a vari- 
able which  we  have  not  yet  been  able  to  change  experimentally  because 
of  fabrication  difficulties  with  lower  Nd  concentration  laser  rods. 

The  approach  we  have  taken  is  straightforward.  Given  the  lamp  emission  spectrum  for  a 
certain  pump  energy,  the  laser  rod  excitation  is  calculated,  and  the  coupled  photon  and  popula- 
tion inversion  rate  equations  are  solved  in  the  low-gain  approximation,  with  and  without  the 
presence  of  a saturable  absorber  in  the  laser  cavity.  Many  forms  of  the  rate  equations  are 
found  in  the  published  literature,  but  we  shall  derive  the  equations  used  in  our  numerical  anal- 
ysis, indicating  the  approximations  used  and  our  nomenclature. 

One  simplification  which  we  use  is  the  assumption  of  a uniform  laser  mode  profile.  This 
probably  corresponds  closely  to  the  actual  situation  of  either  a highly  multimode  output,  where 
the  intensity  of  many  different -order  Hermite-Gaussian  modes  gives  a quasi-uniform  distribu- 
tion across  the  laser  rod,  or  of  the  single  lowest  order  mode  of  an  unstable  resonator,  which 
is  designed  with  the  purpose  of  providing  a nearly  uniform  mode  across  the  gain  medium.  The 
effects  of  having  a given  nonuniform  mode  profile  could  be  included  in  the  calculations,  but  with 
much  more  complicated  and  time-consuming  numerical  integration. 

The  next  approximation  concerns  the  nature  of  the  flash -lamp  excitation.  For  numerical 
purposes,  the  time  dependence  of  the  light  output  is  characterized  by  a Gaussian  pulse.  The 
calculated  shape  for  a simple  LC-flash-lamp  circuit  is  somewhat  different,13  but  the  presence 
of  a saturating  inductance  in  our  series  trigger  transformer  smoothes  out  the  pulse  shape  so 
that  it  does  superficially  resemble  a Gaussian.  In  any  event,  this  should  not  be  a critical  factor 
in  determining  the  integrated  output  energies.  Another  simplification,  which  should  be  fairly 
accurate,  is  the  characterization  of  the  lamp  output  by  a blackbody  spectral  distribution.  With 
the  small-bore  lamps  being  used,  large  peak  current  densities  on  the  order  of  7000  A/ cm2  pro- 
duce a nearly  blackbody  continuum  output,  rather  than  the  narrow  line  emission  found  at  low 
current  densities.  In  principle,  the  temperature  used  in  calculating  the  spectral  distribution 
varies  during  the  pulse,  giving  a distribution  varying  with  time.  To  avoid  this  complication, 
we  assume  that  the  color  temperature  (determined  by  the  electrical  input  energy  and  pulse  dura- 
tion) is  constant,  and  that  the  instantaneous  output  power  is  proportional  to  the  discharge  Gauss-  j 

ian  pulse  amplitude. 
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Fig.  10.  Simplified  energy-level  diagram 
of  Nd^  + indicating  the  levels  used  in  the 
rate  equations.  The  crystal-field  split 
manifolds  are  indicated  as  single  levels, 
and  the  separations  are  not  to  scale. 


2.  Derivation  of  Laser  Rate  Equations 

The  Nd3  + ion  is  modeled  by  the  simplified  energy  level  diagram  in  Fig.  1 0.  The  population 
rate  equation  is 

dN2  ^2  ^2  i 

dt  ‘ “ t “I  r + I B T + /JNo  (21) 

s s a 


where 


N2  = r 3/2  population  density 


4 

Nj  = population  density 


Nq  = ground-state  population  density  “N  = Nd  concentration 
PNq  = excitation  rate  (per  unit  volume) 


t = 4F, 


fluorescence  lifetime 


I = laser  mode  intensity 

Ig  = laser  saturation  intensity  = h v/a  t 

d 4 4 

a = population  weighted  Fj^2  — Ij^2  emission  cross  section 

A 4 4 

a = population  weighted  “*  ^i/z  a^aorPtion  cross  section. 

In  reducing  the  equations  for  the  three  levels  (0),  (1),  and  (2)  to  the  single  necessary  equation 

for  N,,  we  have  used  the  fact  that  levels  (1)  and  (2)  are  thermalized  by  the  phonon  bath,  so  that 
z -E./kT 

Nj  Nq  e (where  Ej/kT  “ 10)  and  that  N2  « NQ,  so  that  the  ground  state  is  only  slightly 

depleted  during  lasing.  All  of  the  levels  shown  are  actually  crystal-field  split  manifolds,  and 


**  *1*  ' ■■  v'1  ’ v"'^  jg. 


«^wWU'. 


the  cross  sections  9 and  o are  the  peak  emission  and  absorption  cross  sections  found  from 
a population -weighted  sum  over  transitions  between  different  levels  in  each  manifold.  Both  N. 
and  0Nn  are  spatially  varying  quantities,  depending  on  the  pump  absorption  parameters  as  de- 
scribed  later.  Defining  the  inversion  density  as  n = N,  - N.(o  /o  ),  Eq.  (21)  becomes 


To  remove  the  spatial  dependence  from  the  inversion  n,  we  integrate  over  the  laser  mode  vol 
ume.  Here  we  use  the  approximation  of  a uniform  laser  mode  intensity  I and  also  neglect  the 
longitudinal  variation  of  I arising  from  spatial  hole-burning  effects.  Defining 


we  have 


0NdV  - V. 


volume  of  laser  mode  in  rod.  We  define  the  pumping  rate  by 


0NdV  - V 


giving 


where  we  introduce  the  parameter  Tq,  the  cavity  lifetime,  for  later  normalization  purposes. 
In  terms  of  other  laser  parameters,  tq  = [2£/c(L  + T)J,  where 


£ = laser  cavity  optical  length 
c = speed  of  light 
L = internal  laser  loss/round  trip 
T = output  mirror  transmission. 

The  next  rate  equation  that  is  needed  is  for  the  saturable -absorber  Q -switch  in  the  laser 
cavity.  In  Fig.  11,  taken  from  earlier  studies  on  mini -laser  rangefinders,14  is  shown  a simpli- 
fied energy-level  diagram  for  the  absorbing  dye  bis-(4-dlmethylaminodithiobenzil)-nickel  or  the 
analogous  ethyl  complex,  abbreviated  as  BDN.1 5 The  important  parameters  for  the  dye  rate 
equation  are  o D>  the  absorption  cross  section  from  the  singlet  ground  state,  and  td>  the  decay 
time  from  the  singlet  excited  state.  There  is  absorption  from  the  excited  singlet  Sj  state  to  the 
second  excited  singlet  S2>  but  since  the  nonradiative  relaxation  back  to  Sj  is  so  fast,  the  popula- 
tion of  St  is  presumed  to  be  unaffected.  The  St  — S2  absorption  with  cross  section  o e will  be 
included  later  as  a loss  term  for  the  laser  photons.  With  these  approximations,  the  simplified 
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Fig.  1 1 . Energy-level  diagram  for  BDN 
dye,  showing  ground  and  excited-state 
cross  sections,  and  o , and  excited- 

state  decay  time,  Tp.  The  absorptions 
indicated  are  at  the  lasing  energy. 
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rate  equations  for  the  ground-state  population  density,  nij,  and  the  first  excited  singlet  density, 


Note  that  nij  + m2  = mQ,  the  dye  molecule  density.  In  terms  of  the  population  difference,  m 
iflj  - m2,  and  the  dye  density  m , m,  = (n-  + m)/2  and  m2  = (m  - m)/2,  giving 


Normalizing  the  variables  with  respect  to  the  dye  density,  saturation  intensity,  and  cavity  life 
time,  Eq.  (27)  becomes 


As  a final  check  of  Eq.  (28),  we  note  that  in  the  state 


in  agreement  with  the  steady-state  populations  calculated  in  Ref.  14. 

The  final  variable  for  which  we  derive  a rate  equation  is  the  laser  mode  Intensity  or,  equiv- 
alently, the  number  of  laser  photons,  q,  within  the  cavity.  If  standing-wave  effects  are  neglected, 
for  a uniform  laser  mode  of  area  Am»  q = (l/l8)  Am£/(cc®T).  Neglecting  spontaneous  emission 
into  the  laser  mode,  the  photon  number  rate  equation  is 
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where 


1 D = Q -switch  length 


A_  = laser  mode  area 
m 


Expressing  q in  terms  of  the  normalized  intensity  gives 
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In  Eq.  (31)  we  choose  to  normalize  the  integrated  inversion  by  a parameter  Jlj,  = Am(L  + T)/2 »B, 
which  would  be  the  threshold  inversion  in  the  absence  of  the  saturable  absorber.  In  terms  of 
this  parameter,  we  also  define  RT  = Tl^/  t,  which  is  the  threshold  pump  rate  without  dye  Iobs. 
With  these  definitions,  Eq.  (25)  becomes 


® (4) = 4 ^ Hi)  4 " 4 + 

and  Eq.  (31)  becomes 
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Next,  we  express  time  in  terms  of  a dimensionless  variable  t/ Tq,  and  define  the  normalized 
variables  by 

)t  3 ji/nT 

1 ■ Vi, 

m = m/mQ 

giving  the  final  dimensionless  rate  equations: 

(gain  inversion)  fi  = A(  —I  j|  -ft  +p]  (34) 
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(intensity) 


(35) 


I = In  — T(1  +|  (1  +«)]  -Im|  (1  -€) 

(dye  population)  m = r[-Iy(m  + 1 ) - in  + 1 ) , (36) 

where  the  dimensionless  parameters  are 
X = tq/t 

p = R/Rt 

P = 2mQiDaIy/(L  + T)  = ratio  of  unsaturated  dye  loss  to  cavitv  loss 
€ = ae/°D 

r = Vtd 

/ B 

y = oD Tj/a  t 

and  the  dot  superscript  denotes  differentiation  with  respect  to  the  normalized  time  variable. 

The  final  quantity  which  is  necessary  for  solution  of  the  rate  equations  is  the  normalized 
pumping  parameter,  p,  which  we  now  calculate  using  the  lamp  emission  and  rod  absorption. 

The  actual  laser  heads  we  have  used  have  a close  coupling  geometry,  with  the  lamp  and  square 
rod  in  close  proximity,  and  both  surrounded  by  an  aluminized  circular  cross  section  reflector. 
This  arrangement  leads  to  a nonuniform  and  complex  absorption  profile,  so  for  ease  of  calcula- 
tions we  have  chosen  to  model  the  arrangement  using  a cylindrical  rod  with  uniform  illumination 
around  the  circumference.  The  intensity  of  the  blackbody  source  at  the  rod  surface  is  assumed 
to  be 


Po(X)  dX  = nL2irh 


c2 

7 


dX 

ehc/XkT_1 


(37) 


where  is  a lamp  radiative  efficiency  factor.  We  will  assume  that  the  flux  in  the  rod  at  each 
wavelength  is  exponentially  decreasing  as  a function  of  the  radius,  so  that  the  volume  integrated 
pumping  rate  is  found  from  integrating  the  spatially  dependent  power  density  divided  by  the  pump 
photon  energy,  and  multiplying  by  the  absorption  coefficient,  (/3N)^  = opr(X)/hv,  giving 


where 


a(X)  = NPP  absorption  coefficient 
rQ  = laser  rod  radius 


r^j  = laser  mode  radius 

For  generality  we  have  allowed  the  rod  and  mode  radii  to  be  different,  so  that  we  may  examine 
the  effects  of  having  incomplete  filling  of  the  rod  volume  by  the  laser  mode.  The  absorption 
coefficient  is  an  experimentally  measured  quantity.  The  electrical  input  energy  to  the  flash 
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lamp  appears  implicitly  in  Eq.  (3  8)  through  the  blackbody  temperature  T.  This  temperature  is 
found  by  equating  the  peak  radiated  power  to  the  integrated  spectral  intensity. 


Plamp  = ffT 


where 


a - Stefan -Boltzmann  constant  = 5.67  x 10"8W/m2K4 
rQ  = lamp  discharge  radius  = laser  rod  radius 

t = lamp  length  = laser  rod  length  , 

and  then  equating  the  integral  of  times  a normalized  Gaussian  time  pulse  to  the  electrical 

energy  into  the  lamp.  Finally,  the  resonant  absorption  term  in  Eq.  (24)  is  subtracted  from 
Eq.  (38)  to  give  the  pumping  rate  R,  which  is  then  normalized  by  R,p.  We  have  also  included 
an  attenuation  factor  of  0.5  to  allow  for  the  effects  of  pump  collector  absorption  by  the  aluminum 
coating  (about  10%)  and  for  obscuration  by  the  lamp  itself  and  imperfect  lamp  imaging  onto  the 
rod.  The  coupled  Eqs.  (34)— (36 ) were  solved  numerically  on  an  IBM  370/168  computer  using  a 
standard  Runge-Kutta  method.  The  solutions  were  carried  out  with  the  normalized  values,  and 
then  scaled  to  the  actual  physical  values  using  the  input  parameters  such  as  pump  energy,  cavity 
loss,  and  saturable-absorber  loss.  Normal  mode  (non-Q -switched)  laser  behavior  was  found 
quite  simply  by  setting  the  parameter  /3  = 0. 

Up  to  this  point  we  have  treated  the  Nd  concentration  as  being  a fixed  parameter,  N = 3.9  x 
10Z*  cm  3 for  pure  NdP^O^.  However  it  is  possible  to  replace  Nd  by  the  optically  inert  rare- 
earth  La  to  obtain  Nd^a^^PjOj^,  so  that  one  has  control  over  some  of  the  fluorescence  prop- 
erties of  the  laser  material.  The  primary  effects  of  varying  the  fractional  Nd  concentration, 
x,  are  threefold;  (1)  The  absorption  coefficients,  approximately  proportional  to  x,  will  change 
and  give  different  spatial  inversion  and  gain  densities  across  the  rod  cross  section.  (2)  The 

4 

thermal  population  of  the  terminal  Ij  j/2  manif°ld  will  vary  with  x,  changing  the  resonant  self- 
absorption loss  and  laser  threshold.  (3)  The  nonradiative  concentration  quenching  of  the  Nd3  + 
fluorescence  decreases  at  lower  Nd  concentration,  thus  increasing  the  fluorescent  lifetime,  t, 
of  the  upper  4Fj/2  laser  manifold  as  x decreases.  All  of  these  effects  can  be  readily  incor- 
porated in  the  rate  equations  by  making  the  following  substitutions: 


N — Nx 

T T0/<1  + qx)  ' 

where  rQ  = 350  psec  and  qQ  =■  1.69,  as  determined  from  our  own  experimental  data.  Other 
numerical  values  which  were  used  in  the  calculations  are  as  follows: 

£ = laser  cavity  length  = 10  cm 

rQ  = rod  and  lamp  radii  = 0.075  cm 

A -VkT  -24  2 

a e = effective  resonant  absorption  cross  section  = 6x10  cm 

d —1 Q 2 

a = effective  emission  cross  section  = 1.5X10  cm 
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Figs.  12-15.  Gain  as  a function  of  time  for  1-cm-long  NPP  laser  rod  with  flash-lamp 
pump  having  the  indicated  relative  time  dependence.  Each  gain  curve  is  for  a different 
fractional  Nd  concentration,  ranging  from  x = 1.0  to  x = 0.25. 
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Ojj  = ground-state  dye  absorption  cross  section  = 1.2  x 10-16  cm2 

ae  = excited-state  dye  absorption  cross  section  = 2.85  x to'1 7 cm2 

-9 

Tjj  = dye  decay  time  = 1.87  x 10  sec 

The  input  variables  for  each  calculation  were  the  excitation  energy.  Nd  concentration,  cavity 
loss  (assumed  to  be  entirely  from  the  output  coupler),  unsaturated  dye  loss,  and  mode  area. 

3.  Absorption  and  Gain  Characteristics 

Before  giving  the  calculated  laser  results,  we  examine  some  of  the  absorption  and  gain  char- 
acteristics. In  Figs.  12  to  15  is  shown  the  gain  2gt  as  a function  of  time  for  1 J pump  energy  in 
lasers  with  different  Nd  concentration  x = 1.0,  0.75,  0.50.  0.25.  The  gain  parameter  plotted  on 
the  vertical  axis  is  the  round-trip  gain  of  the  uniform  laser  mode  through  the  1-cm-long  laser 
rod,  and  is  found  by  integrating  the  spatially  nonuniform  inversion  over  the  rod  cross  section. 

As  expected,  the  gain  decreases  with  lower  absorption  and  the  peak  of  the  gain  shifts  slightly  to 
longer  times  as  the  fluorescence  lifetime  increases.  However,  for  the  relatively  short  pump 
pulse  being  used  (with  normalized  time  dependence  given  by  the  dashed  curves),  the  pump  is 
absorbed  in  a time  small  compared  with  decay  or  quenching  times  for  any  concentration.  There- 
fore. even  though  the  fluorescence  decay  characteristics  for  times  greater  than  100  psec  can 
be  seen  to  differ  with  different  values  of  x,  the  gains  for  shorter  intervals  do  not  depend  on  t 
directly,  but  rather  on  the  pump  absorption  strength  and  to  some  extent,  on  the  resonant  absorp- 
tion loss.  Clearly  this  will  not  be  the  case  for  C'W  lasers,  where  r plays  a direct  and  signifi- 
cant role  in  determining  the  steady-state  inversion  and  gain.  Another  interesting  point  seen  in 
these  curves  is  that  the  peak  of  the  gain  occurs  30  to  40  psec  after  the  peak  of  the  pump  pulse, 
so  that  an  optimally  designed  Q-switched  laser  should  hold  off  lasing  until  this  maximum  energy 
storage  condition  is  reached  in  order  to  extract  the  most  energy  in  a single  short  pulse. 

The  spatial  dependence  of  the  pumping  found  from  the  absorption  spectrum  of  NPP  and  the 
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blackbody  lamp  emission  is  shown  in  Fig.  16,  for  a 1 -J  pump  energy  and  the  four  previous  Nd 
concentrations.  Because  of  the  large  peak  absorption  coefficients  in  NPP,  the  strongest  con- 
centration dependence  on  the  inversion  appears  near  the  surface  of  the  laser  rod.  In  the  interior, 
the  inversion  is  more  insensitive  to  concentration  because  the  weaker  peak  absorptions  of  the 
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dilute  (Nd,  La)PP  give  absorption  strengths  comparable  with  those  arising  from  the  wings  of 
the  absorption  distribution  in  NPP.  For  the  uniform  laser  mode  the  differential  inversion, 
m(r),  is  plotted  in  Fig.  1 7,  and  shows  the  relative  contribution  to  the  mode  gain  from  the  inver- 
sion at  a distance  r from  the  rod  center. 

4.  Normal- Mode  Laser  Output 

Next,  a typical  laser  output  pulse  is  shown  in  Fig.  18,  in  this  case  with  x = 1.0  and  1 J pump 
energy  into  the  lamp.  We  have  taken  the  output  mirror  transmission  (representing  the  entire 
cavity  loss)  equal  to  0.2  and  have  not  included  any  dye  loss.  The  oscillatory  output  comes  from 
the  relaxation  oscillations  induced  by  the  rapid  population  inversion  and  overshoot  from  the 
equilibrium  threshold  level.  Such  oscillations  are  commonly  observed  in  solid-state  lasers, 
particularly  those  having  long  fluorescent  lifetimes.  In  our  case,  the  regularity  of  the  oscilla- 
tion comes  from  our  assumption  of  a single-mode  output.  In  actual  practice,  multi -longitudinal 
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Fig.  16.  Inversion  density  n(r)  as  a function 
of  radial  coordinate  in  a 1.5-mm-dia.  (Nd, 
La)  PP  rod  with  1 J pump  energy  for  the  dif- 
ferent fractional  Nd  concentrations  indicated. 


Fig.  17.  Differential  inversion  density 
rn(r)  as  a function  of  radial  coordinate 
in  a 1.5-mm-dia.  (Nd,  La)  PP  rod  with 
1 .)  pump  energy  for  the  different  frac- 
tional Nd  concentrations  indicated. 
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Fig.  18.  Calculated  normalized  output  intensity  i/l 
as  a function  of  time  for  a 1.5-mm-dia.,  l-cm-lon| 
NPP  rod  with  1 J energy  input  and  20  percent  output 
mirror  transmission.  The  integrated  absolute  out- 
put energy  is  24  mj. 
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Fig.  19.  Calculated  output  energy  as  a function  of  input 
energy  (dotted  curve,  top  scale)  for  the  NPP  laser  with 
complete  mode  filling  of  the  rod  (r^j  = rQ).  Output  en- 
ergy as  a function  of  Nd  concentration  x with  1 J energy 
input,  and  r^j  = rQ  (solid  curve,  bottom  scale)  and  r^j  = 

3/4  rQ  (dashed  curve,  bottom  scale). 

and  transverse  modes  cause  a highly  irregular  spiking.  In  any  case,  we  are  not  concerned  so 
much  with  the  detailed  time  behavior  of  the  normal-mode  output,  but  rather  its  time  integral, 
which  gives  the  output  pulse  energy. . In  this  case,  we  have  an  output  energy  of  .24  mj,  with  1 J 
input. 

Curves  and  data  for  normal -mode  operation  similar  to  Fig.  18  were  calculated  varying  sev- 
eral parameters  and  these  are  summarized  in  Fig.  19.  First,  as  a check  on  the  validity  of  the 
calculations,  we  computed  a curve  of  output  vs  input  energies  (dotted  curve,  top  horizontal  axis). 
Both  the  threshold  and  peak  output  energy  are  in  good  qualitative  agreement  with  experimental 
results  to  be  discussed  in  the  next  section.  The  slight  deviation  from  linearity  probably  stems 
from  the  conversion  of  input  energy  to  blackbody  lamp  temperature,  with  a consequent  lower- 
ing of  the  emission  peak  from  0.54  |tm  at  ~0.2  J input,  to  0.36  jim  at  1 J input.  Such  a shift  in 
the  emission  spectrum  away  from  the  Nd  absorption  regions  would  cause  a nonlinear  decrease 
in  absorption  efficiency.  Next,  to  compute  the  optimum  Nd  concentration  for  maximum  energy 
output  we  calculated  EQut  as  a function  of  concentration  x (lower  horizontal  axis)  for  a fixed  in- 
put energy  of  1 J (solid  curve).  As  might  be  expected  from  the  earlier  gain  calculations,  the 
maximum  concentration  x = 1.0  gives  the  best  absorption  and  the  most  output  energy,  with  the 
smaller  effects  of  larger  concentration  quenching  and  resonant  loss  not  being  significant.  The 
calculations  for  the  very  lowest  concentrations  (x  & 0.3)  are  probably  not  as  accurate  as  the 
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high-concentration  values  because  our  laser  model  did  not  include  reabsorption  of  pump  radia- 
tion partially  transmitted  through  the  laser  rod  and  re-reflected  by  the  pump  collector  back  Into 
the  rod. 

In  view  of  the  results  of  Figs.  16  and  17  which  indicate  nonuniform  absorption,  we  have  cal- 
culated the  effects  of  having  a laser  mode  which  does  not  completely  fill  the  rod,  as  might  be 
caused  by  edge  diffraction  loss,  for  example.  The  dashed  curve  in  Fig.  19  shows  the  calculated 
result  for  r^j  = 0.75  r , in  this  case  a 1.13-mm-dia.  mode  in  a 1.5-mm-dia.  rod.  The  maxi- 
mum output  decreases  from  the  previously  mentioned  24  to  8 mj.  Even  though  the  concentra- 
tion dependence  is  significantly  reduced,  the  largest  output  energy  is  still  obtained  for  x = 1. 

5.  Q-Switched  Laser  Output 

When  the  saturable  absorber  is  inserted  into  the  cavity,  lasing  oscillation  is  suppressed 
until  the  gain,  as  shown  in  Fig.  12  for  example,  exceeds  the  combined  cavity  and  unsaturated 
dye  loss.  The  laser  oscillation  then  builds  up  rapidly  in  a time  on  the  order  of  Tq,  and  as  it 
does,  the  laser  mode  intensity  bleaches  the  dye,  lowering  the  total  loss,  and  increases  the  mode 
buildup  still  further.  This  process  continues  until  the  population  inversion  is  depleted  and  the 
mode  decays,  resulting  in  a single  short  output  pulse  if  the  pump  and  loss  parameters  are  prop- 
erly chosen.  If  the  dye  loss  is  too  large,  no  lasing  will  occur,  and  if  it  is  too  small,  multiple 
output  pulses  may  result. 


Fig.  20.  Calculated  Q-switched  output 
intensity  as  a function  of  time  for  NPP 
rod  with  i J input.  The  output  coupling 
is  0.2,  the  unsaturated  dye  absorption 
is  0.8,  and  the  integrated  energy  is 
10.2  mj. 


A typical  Q-switched  output  pulse  for  1 J input  and  x = 1 is  shown  in  Fig.  20.  The  output 
mirror  transmission  is  0.2  and  the  unsaturated  dye  loss  is  0.8,  near  the  maximum  allowed  for 
1 J input.  The  time  scale  is  shown  greatly  expanded  compared  with  Fig.  18,  and  only  a small 
portion  containing  the  single  output  pulse  ~75  fiBec  after  the  start  of  the  pump  is  shown. 

All  parameters  were  held  fixed  (Em  = 1 J,  x = 1.0)  except  for  mirror  transmission  and  un- 
saturated dye  loss,  and  the  resulting  Q-switched  output  behavior  Is  summarized  In  Fig.  21.  First, 
with  the  output  mirror  transmission  held  at  T = 0.2,  the  Q-swltch  loss  was  varied  from  0.45  to 
0.8,  approximately  the  range  between  threshold  (T  £ 0.85)  and  double -pulsing  (T  £ 0.4).  The 
output  energy  (open  circles,  dashed  curve,  left  scale)  Increased  monotonically  from  7.8  to 
10.2  mj,  and  the  pulse  width  (full-width,  half-maximum)  decreased  monotonically  from  9.0  to 
5.6  nsec  (open  circles,  solid  curve,  right  scale).  These  results  showed  that  at  1 J input,  the 
maximum  total  loss  was  1.0,  so  this  parameter  was  fixed  and  both  the  Q-switch  loss  and  mirror 
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Fig.  21.  Calculated  Q-switched  energy 
output  (left  scale,  dashed  curves)  and 
full-width,  half-maximum  pulse  dura- 
tion (right  scale,  solid  curves)  as  a func- 
tion of  Q-switch  loss.  The  curves  con- 
necting the  circles  refer  to  fixed  output 
mirror  coupling  of  0.2,  and  the  curves 
connecting  the  squares  refer  to  variable 
coupling,  maintaining  the  total  coupling 
plus  Q-switch  loss  at  1.0. 


transmission  were  varied,  but  maintaining  the  sum  constant  at  1.0.  This  procedure  showed  that 
the  maximum  theoretical  energy  extraction  occurred  for  the  optimum  combination  of  T = 0.4  and 
Q-loss  = 0.6,  with  an  energy  output  of  12.1  mJ  and  a pulse  width  of  6.5  nsec.  The  maximum  in 
the  power  dependence  is  a weak  one,  but  it  would  be  better  to  choose  that  operating  point  pri- 
marily for  the  most  stable  laser  output,  midway  between  the  threshold  and  double -pulsing  con- 
ditions. 


B.  EXPERIMENTAL  STUDIES 

1.  Overview 

The  experimental  studies  of  NPP  lasers  that  we  have  undertaken  are  a continuation  of  our 
original  work  using  flash-lamp  excitation,  which  was  previously  reported.4  The  main  effort  was 
to  increase  the  output  energy  and  efficiency  of  the  laser,  to  study  its  Q-switched  behavior  and 
maximize  the  short-pulse  output  energy,  and  in  so  far  as  possible,  provide  experimental  evi- 
dence to  give  confidence  in  our  numerical  simulation,  described  above. 

One  of  the  laser  parameters  which  we  had  varied  in  the  computer  calculations  and  which  we 
hoped  to  study  experimentally  as  well  was  the  Nd  concentration  in  the  laser  rod.  Although  we 
expect  that  pure  NPP  will  have  the  best  absorption  efficiency,  other  factors  such  as  the  laser 
mode  nonuniformity  and  the  details  of  the  absorption  profile  in  the  square  rod  make  the  laser 
output  calculations  somewhat  imprecise.  Clearly,  it  would  be  desirable  to  have  rods  of  (Nd,  La) 
PP  to  study  such  effects  and  compare  with  the  calculations  of  Fig.  19.  Unfortunately,  as  La 
replaces  Nd,  large  high-quality  pentaphosphate  crystals  become  more  difficult  to  grow.  The 
main  problem  seems  to  arise  from  strains  in  the  mixed  crystals  which  make  it  nearly  impos- 
sible, at  the  present  time,  to  cut  and  polish  laser  rods  without  introducing  cleavage  cracks. 
Although  further  studies  of  the  growth  process  may  improve  this  situation,  because  of  the  un- 
availability of  the  dilute  rods,  all  of  the  experimental  results  described  in  this  report  were 
obtained  with  pure  NdP5014  laser  rods. 
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2.  Normal-Mode,  Stable  Resonator  Operation 

In  this  section,  we  describe  the  experimental  results  using  a stable  confocal  resonator  (with 
10  cm  mirror  separation)  and  normal-mode  (non-Q-switched)  operation. 

The  laser  parameters  which  were  investigated  under  these  conditions  were  the  effects  of 
rod  size,  end  facet  coatings,  and  pumping  uniformity.  The  basic  laser  system  was  similar  to 
that  described  in  Ref.  4.  The  laser  cavity  consisted  of  dielectric -coated  mirrors  with  10-cm 
radii  of  curvature.  One  mirror  had  high  reflectivity  (>99.9%)  at  1.05  pm,  and  different  end  mir- 
rors with  reflectivities  from  97.5%  to  80%  were  used.  The  laser  head  consisted  of  an  aluminum 
body  which  was  drilled  out  to  accommodate  a 6-mm-i.d.  glass  tube,  whose  inner  surface  was 
aluminized  by  evaporation  from  an  aluminum -wetted  tungsten  filament.  One  edge  of  the  reflect- 
ing cylinder  was  slotted  parallel  to  its  axis  to  allow  insertion  of  a metal  tab  to  which  was  ce- 
mented the  edge  of  the  laser  rod.  The  side  of  the  aluminum  body  was  also  slotted  next  to  the 
reflector  slot,  and  the  mounting  tab  was  clamped  by  the  aluminum  body. 

The  flash  lamps  had  hard  borosilicate  glass  envelopes,  with  a 1.5-mm  i.d.,  and  a high- 
pressure  fill  of  Xe  gas.  The  electrode  separation  was  varied  to  match  the  length  of  the  laser 
rod,  but  in  most  cases  was  1 cm.  These  lamps  are  an  experimental  type,  supplied  to  us  by  the 
GTE  Sylvan ia  Lighting  Products  Group.  Because  the  fabrication  procedure  was  newly  developed, 
no  precise  control  or  knowledge  of  the  fill  pressure  was  available,  but  it  was  estimated  to  be 
approximately  1500  Torr. 

The  lamp  was  inserted  in  the  pump  reflector  assembly  parallel  with  and  immediately  adja- 
cent to.  but  not  contacting,  the  laser  rod.  It  was  mounted  by  cementing  the  glass  envelope  to 
plastic  end  pieces  attached  to  the  ends  of  the  aluminum  reflector  body.  Bayonet-type  connector 
pins  were  soldered  via  flexible  copper  sheet  leads  to  the  lamp  electrodes.  The  lamp  was  fired 
through  a standard  LC  discharge  network.  We  used  a iO-uF  (2  kV  rated)  storage  capacitor  with 
a 60-pH  coil  inductance.  The  lamp  trigger  voltage  was  supplied  through  a series  trigger  trans- 
former unit  (EGG  type  TM-12A),  which  added  a small,  saturating  inductance  to  the  circuit.  The 
current  waveform  was  a smooth,  bell -shaped  curve  with  about  a 50-nsec  FWHM  duration,  and 
peak  current  (at  1 J)  of  about  120  A.  The  repetition  rate  was  set  with  a low-frequency  function 
generator,  which  in  turn  triggered  smother  pulse  generator  with  pulse-delay  capability,  shorter 
pulses,  smd  higher  output  voltages.  This  narrow-pulse  output  triggered  the  EGG  high-voltage 
trigger  module,  discharging  the  10-uF  capacitor  which  had  been  charged  with  a standard  variable - 
voltage  DC  supply  through  a 30-kfl  current -limiting  resistor. 

The  laser  output  pulses  were  measured  with  a Si  photodiode  terminated  by  50  f)  (for  fast  re- 
sponse) at  the  oscilloscope  input.  Integrated  pulse  energies  were  measured  by  a Scientech  362 
thermopile  meter,  and  generally  obtained  from  an  average  of  sev>  al  laser  shots,  i.e.,  measur- 
ing the  total  energy  of  about  ten  successive  shots  and  dividing  by  the  number  of  shots.  For  ex- 
amination of  the  laser  transverse  mode  distribution,  a flip-in  mirror  diverted  the  laser  beam 
onto  an  opaque  screen,  at  the  focus  of  an  infrared-sensitive  TV  camera  whose  picture  was  dis- 
played on  a TV  monitor. 

First,  we  shall  review  the  results  obtained  with  rods  of  different  sizes.  All  rods  described 
below  had  the  "c*  axes  parallel  to  the  lasing  direction,  and  had  polished  end  faces.  The  most 
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significant  improvement  over  the  initial  results  obtained  with  the  same  size  rod  came  from 
four  factors:  (1)  adjusting  the  end  mirrors  for  multimode  output  in  order  to  fill  the  rod  volume, 

(2)  using  high-pressure  Xe  lamps  for  more  efficient  operation,  (3)  using  a slightly  more  closely 
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Fig.  22.  Output  vs  input  energy  for  a Fig.  23.  Output  vs  input  energy  for  a 

1-  x 1-  x 7 -mm  NPP  rod  in  a 10-cm  i-  x 2-  x 10-mm  NPP  rod  in  a 10-cm 

confocal  cavity  with  9.5%  and  20%  trans-  confocal  cavity  with  9.5%  and  20% 

mitting  output  mirrors.  transmitting  output  mirrors. 

coupled  pump  collector,  and  (4)  lowering  the.  repetition  rate  to  12  pulses  per  minute  (0.2  Hz)  to 
eliminate  rod  heating  effects  and  ensure  full  charging  of  the  lamp  capacitor.  With  these  improve 
ments,  we  have  obtained  laser  outputs  shown  in  Fig.  22,  using  a 1-  x 1-  x 7 -mm  rod  and  two  dif- 
ferent output  mirrors  of  9.5%  and  20%  transmission.  In  this  case,  the  slope  efficiencies,  ~1  per- 
cent, are  nearly  the  same,  implying  that  the  internal  laser  losses  are  much  less  than  10  percent. 

Figure  23  shows  the  results  using  a 1-  x 2-  x 10-mm  rod,  with  the  broad  2-  x 10-mm  face 
adjacent  to  the  flash  lamp.  The  maximum  output  energy  at  1 J input  has  increased  nearly  40  per- 
cent to  11  mj,  and  the  slope  efficiency  for  the  20%  T mirror  has  increased  to  ~1.4  percent. 
Again  the  most  significant  difference  between  the  results  with  the  two  output  couplers  is  the  in- 
creased threshold  with  increasing  mirror  transmission,  although  here  there  is  also  a slight  but 
noticeable  increase  in  slope  efficiency  as  well. 

In  order  to  obtain  more  energy  from  the  NPP  laser,  the  next  step  was  to  increase  the  length 
of  both  the  rod  and  lamp.  This  was  done  primarily  in  order  to  drive  the  lamp  to  higher  levels 
without  shortening  its  life.  The  lamp  explosion  energy  for  the  glass,  high-pressure  lamps  is 
not  known,  but  the  value  for  quartz,  low-pressure  (450  Torr)  lamps  is  estimated  from  various 
manufacturers'  literature1^1  to  be  on  the  order  of  20  J/cm  for  our  present  pulse  width.  The 


Fig.  24.  Output  vs  input  energy  for  a 
2-  x 2-  x 20-mm  NPP  rod  in  a 10-cm 
confocal  cavity  with  a 20%  transmitting 
output  mirror. 
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Fig.  25.  Output  vs  input  energy  for  a 
2-  x 2-  x 10-mm  NPP  rod  with  anti- 
reflection coatings  in  a 10-cm  confo- 
cal cavity  with  9.5%  and  20%  trans- 
mitting output  mirrors. 


Fig.  26.  Output  vs  input  energy  for  a 2-  x 
2-  x 10-mm  NPP  rod  with  one  end  having 
a high-reflection  coating  and  the  other  an 
anti-reflection  coating.  The  semiconfocal 
cavity  was  5 cm  long  with  9.5%  and  20% 
transmitting  output  mirrors. 
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lamp  lifetime  (number  of  flashes)  should  increase  nearly  exponentially  with  decreasing  operat- 
ing energy,  expressed  as  a percentage  of  the  explosion  energy.  In  our  case,  the  factors  of  a 
higher  pressure  fill  and  a less  durable  envelope  would  probably  decrease  the  explosion  energy, 

possibly  to  ~10  J/cm.  Therefore,  we  estimate  that  for  a lamp  input  of  4 J,  the  life  of  a longer 
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lamp  would  be  extended  from  ~10  to  10  flashes  to  >10  flashes.  Another  benefit  would  be  to 
lower  the  peak  current  density  and  color  temperature,  shifting  the  blackbody  emission  peak 
towards  the  Nd3+  absorption  bands.  Using  a larger  rod  cross  section  of  2 x 2 mm  and  the  pre- 
vious lamp  diameter,  but  increasing  both  rod  and  lamp  to  2 cm  length  gives  the  results  shown 
in  Fig.  24  with  a 20%  T output  mirror.  In  this  case,  the  longer  lamp  led  to  unreliable  firing 
below  ~350  V,  so  that  the  laser  threshold  and  lowest  level  outputs  were  not  measured.  Extrap- 
olating from  the  data,  the  estimated  threshold  is  ~400  mj.  The  output  remains  linear  up  to 
input  energies  near  3.4  J,  with  a slope  efficiency  of  2.5  percent.  In  this  regime,  we  obtain  a 
maximum  energy  of  72  mj  with  3.4  J input,  for  a 2.1 -percent  total  efficiency.  Beyond  this  en- 
ergy, the  slope  efficiency  decreases  markedly  to  1.3  percent,  and  the  maximum  output  at  4 J is 
80  mj,  whereas  the  output  with  continued  2.5-percent  slope  efficiency  would  be  88  mJ.  The  most 
probable  cause  for  the  decrease  in  efficiency  is  related  to  the  onset  of  oscillating  modes  not  re- 
quiring the  80% -reflecting  output  mirror.  We  observed  that  with  removal  of  the  output  mirror, 
lasing  occurred  at  a threshold  near  3.2  J.  This  means  that  sufficient  gain  had  built  up  to  allow 
oscillation  with  the  high-reflecting  end  mirror  and  the  5%  reflectivity  of  the  uncoated  NPP  rod 
facet.  The  near-field  output  of  this  mode  was  a narrow  stripe  with  orientation  parallel  to  the 
pumped  face  of  the  laser  rod,  indicating  that  sufficient  gain  for  this  mode  occurred  only  near 
the  pumped  surface.  Using  the  high-gain  approximation,  with  the  threshold  condition  exp  (2gl)  = 
l/R,  gives  the  gain  coefficient  g ~ 0.75  cm  . 

In  all  of  the  results  given  so  far,  the  ends  have  been  polished  but  no*  coated.  In  order  to  see 
if  these  uncoated  faces  were  causing  any  significant  losses  in  the  laser,  we  have  measured  the 
laser  output  using  a rod  with  dielectric  anti-reflection  coatings  on  the  ends.  A 1-  x 2-  x 10-mm 
AR-coated  rod  was  tested  in  the  confocal  cavity  with  the  20%  T output  mirror.  The  threshold 
was  250  mJ,  and  the  output  at  1 J input  was  10.1  mJ.  These  results  are  nearly  identical  with 
those  found  using  the  uncoated  1-  x 2-  x 10-mm  rod  (Fig.  23)  with  threshold  of  250  mJ  and  output 
of  10.9  mJ. 

The  results  with  a 2-  x 2-  x 10-mm  AR-coated  rod  are  shown  in  Fig.  25.  Both  the  overall 
and  slope  efficiencies  are  somewhat  larger  than  in  the  case  of  the  1-  x 1-  x 10-mm  rods.  The 
9.5%  T and  20%  T mirrors  give  slope  efficiencies  of  1.9  and  2.2  percent,  respectively.  As  seen 
in  Fig.  25,  the  maximum  output  energy  is  34.5  mj  with  2 J input  energy.  Over  some  of  the  range 
of  input  energy,  there  seems  to  be  a systematic  deviation  from  the  linear  output  vs  input  char- 
acteristic, but  it  is  not  known  whether  this  is  associated  with  the  lamp  behavior  or  the  laser-rod 
absorption. 

Next,  a similar  2-  x 2-  x 10-mm  rod  was  tested  with  one  face  AR-coated  as  before,  and  the 
other  face  coated  with  a high -reflectance  dielectric  coating  to  compare  such  an  end  reflector  with 


the  external,  10-cm-radius  HR  mirror.  The  external  HR  end  mirror  was  removed,  leaving  a 
semiconfocal  cavity  with  ~ 5 cm  separation  between  the  HR-coated  rod  face  and  the  external  out- 
put mirror.  From  symmetry  considerations,  the  modes  of  this  semiconfocal  cavity  should  not 
differ  significantly  from  those  of  the  longer  confocal  cavity.  The  output  vs  input  energy  data  for 
this  laser  are  shown  in  Fig.  26,  and  differ  little  from  those  of  the  2-  x 2-  x 10-mm  AR  rod  with 
two  external  reflectors  (Fig.  25). 
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Fig.  27.  Back-illuminated  view  of  1-  x 2-  x 10 -mm  NPP  laser  rod 
with  mounting  tab  below  and  lamp  above.  The  video  analyzer  signal 
from  top  to  bottom  along  the  white  scan  line  is  shown  below  the  laser 
rod  view  with  signal  intensity  proportional  to  brightness. 


A comparison  of  the  experimental  results  with  the  previous  normal -mode  calculations  is 
best  made  using  the  data  of  Fig.  23  for  the  1 - x 2-  x 10-mm  rod.  This  area  of  2 mm  is  close 
to  the  cylindrical  area  of  1.77  mm2  for  the  1.5-mm-dia.  rod  assumed  in  the  numerical  analysis. 
With  the  20%  T output  mirror,  the  experimental  and  theoretical  thresholds  are  250  and  160  mj, 
respectively,  with  output  energies  of  11  and  24  mJ  at  1 J input.  Considering  the  ab  iniUo  nature 
of  the  calculation,  the  agreement  to  within  a factor  of  2 is  quite  good.  The  uncertainty  in  the 
gain  distribution  within  the  rectangular  rod,  and  in  the  exact  mode  irofile  can  easily  account  for 
this  difference,  as  demonstrated  by  the  reduction  of  the  calculated  output  from  24  to  8 mj  with 
a mode  three-quarters  the  original  diameter.  Similarly,  an  overall  reduction  in  the  theoretical 
effective  pumping  efficiency  by  a factor  of  2 would  raise  the  threshold  to  ~300  mJ  and  lower  the 
output  to  12  mj,  giving  excellent  agreement  with  the  experimental  data.  From  these  results, 
we  conclude  that  the  numerical  simulation  is  an  accurate  semi -quantitative  model  which  should 
be  capable  of  predicting  the  concentration  dependence  of  (Nd,  La)  PP  laser  rods,  and  of  modeling 
the  Q -switched  behavior,  described  later. 

Measurements  were  performed  to  determine  the  degree  of  pumping  nonuniformity  in  the  laser 
rod.  The  laser  cavity  mirrors  were  removed,  the  head  rotated  90”  about  the  rod  axis,  and  the 
front  surface  of  the  laser  rod  was  imaged  directly  onto  the  infrared-sensitive  vidicon  tube.  A 
CVI  321  Video  Analyzer  was  inserted  between  the  camera  control  unit  and  the  monitor,  and  pro- 
vided a signal  proportional  to  the  picture  intensity  along  a scan  line  perpendicular  to  the  monitor 
raster  scan.  Figure  27  shows  a photograph  of  monitor  screen  on  which  is  seen  the  face  of  the 
1-  x 2 -mm  laser  rod,  with  the  mounting  tab  below  it,  and  a circular  edge  of  the  pump  lamp  above 
it.  This  picture  was  obtained  with  tungsten  lamp  illumination  of  a card  placed  behind  the  laser 
rod.  The  vertical  white  line  shows  the  scan  axis  for  the  video  analyzer.  Also  shown  in  Fig.  27 
is  the  oscilloscope  display  of  the  analyzer  voltage  along  the  scan  line  and  is  provided  for  cali- 
bration purposes.  The  right  side  of  the  oscilloscope  trace  shows  the  sharp  drop  in  intensity  be- 
tween the  laser  rod  and  dark  mounting  tab  (bottom  part  of  scan),  and  the  left  side  of  the  scan 
shows  a narrow  dip  corresponding  to  the  top  edge  of  the  rod  adjacent  to  the  lamp.  After  obtain- 
ing these  photographs,  the  back  illumination  was  removed,  and  with  attenuators  and  visible- 
blocking,  infrared-transmitting  filters  (Corning  2-64,  7-56)  inserted,  -the  Nd3+  fluorescence 
emission  during  the  flash-lamp  excitation  was  observed.  Figure  28  (top)  shows  a photograph  of 
the  monitor  displaying  the  fluorescence  across  the  rod,  and  two  scans  (bottom)  of  that  fluores- 
cence pattern  in  the  vertical  direction  with  the  right  photograph  showing  a center  scan,  and  the 
left  photograph  showing  a scan  near  the  left-hand  edge  of  the  rod.  As  is  evident  from  Fig.  28, 
there  is  considerable  nonuniformity  of  the  pump  absorption,  particularly  along  the  lamp-rod 
dimension  where  the  absorption  decreases  by  a factor  of  2 across  the  i-mm  rod  thickness. 

(The  multiple  scan  traces  in  Fig.  28  are  present  because  of  the  lack  of  synchronization  between 
lamp  firing  and  monitor  scan,  which  necessitated  photographing  several  analyzer  scans.) 

An  attempt  was  made  to  improve  the  pumping  uniformity  of  the  NPP  rod  by  constructing  a 
dual-lamp  (DL)  laser  head.  A cross  section  comparing  the  single  and  DL  configurations  is 
shown  in  Fig.  29.  As  seen  from  the  sectional  view  of  the  DL  collector,  the  degree  of  coupling 
for  each  lamp  was  less  than  that  in  the  single-lamp  collector.  Although  the  copper  rod  mounting 
tabs  for  the  DL  laser  were  polished  and  aluminized,  we  expect  that  light  from  the  lamps  entering 
the  corners  between  the  tabs  and  the  glass  sleeve  was  not  efficiently  reflected  to  the  laser  rod. 

The  lamps  in  this  unit  were  connected  in  series  to  insure  full  discharge  current  through  both 
lamps.  The  procedure  which  we  followed  to  examine  the  effects  of  pumping  uniformity  was  to 
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Fig.  28.  Fluoresence  emission  from  end  face  of  NPP  laser  rod  (top)  and  video  analyzer 
scans  along  vertical  lines  (bottom).  The  left  bottom  trace  is  a scan  near  the  left  edge 
of  the  rod  view,  and  the  right  bottom  trace  shows  a scan  at  the  rod  center. 
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compare  the  laser  output  with  each  lamp  fired  individually  with  that  when  both  lamps  were  fired, 
maintaining  the  same  energy  input.  Firing  both  lamps  gave  an  output  of  19.2  mJ  at  3 J input  and 
12.5  mJ  at  2 J input.  The  threshold  was  below  the  minimum  firing  level  for  both  lamps  in  series, 
470  V or  1.1  J.  Using  the  2-J  level  for  reference,  one  lamp  gave  an  output  of  7.8  mJ,  and  the 
other  gave  2.3  mJ.  Repeating  the  dual-lamp  measurement  at  1.6  J input  then  showed  a decreased 
output  of  4.9  mJ  from  the  previous  9.1  mJ.  Disassembly  of  the  laser  head  showed  that  the  rod 
had  loosened  from  one  of  the  side  tabs,  possibly  affecting  the  alignment  during  the  test.  Also, 
slight  crazing  on  the  lamp  bore  was  observed,  although  probably  not  enough  to  significantly  alter 
the  light  output.  It  is  possible  that  the  second  lamp  is  indeed  poorer  than  the  first,  but  we  have 
not  had  the  time  to  characterize  all  the  lamps  used,  or  to  repair  the  head  assembly  and  repeat 
the  measurements. 

Another  modification  of  the  single -lamp  laser  head  which  we  investigated  was  the  replace- 
ment of  the  aluminized  reflector  sleeve  by  a coating  of  high-diffuse-reflectance  BaS04  paint 
(Kodak  6080)  directly  on  the  machined  and  etched  aluminum  laser  head  assembly.  Using  the 
2-  x 2-  x 10-mm  AR-coated  laser  rod,  we  obtained  25  mJ  output  at  2 J input,  compared  with 
the  earlier  result  of  34.5  mJ.  A different  lamp  was  used  in  this  test  with  a significantly  higher 
discharge  voltage,  so  the  lamp  may  have  contributed  to  the  lower  output.  In  view  of  the  fairly 
good  laser  output  and  the  simplicity  of  laser  head  fabrication,  this  spray-painted  coating  may 
be  a practical  alternative  to  the  coated  glass  sleeve  reflector. 

In  order  to  demonstrate  the  miniaturization  possible  with  the  flash-lamp-excited  NPP  laser, 
we  have  constructed  a self-contained,  hand-held  laser  unit.  For  increased  ruggedness  and  min- 
iaturization, we  have  replaced  the  hollow  pump  collector  by  a solid  integral  assembly.  This 
consists  of  a 1-cm-long  x 6-mm-dia.  sapphire  cylinder,  with  a round  hole  to  accommodate  the 
2.5-mm-dia.  Xe  flash  lamp  and  a 1-mm-square  slot  to  hold  the  NPP  laser  rod,  both  of  which 
are  mounted  using  polystyrene  cement.  The  outer  surfaces  of  the  sapphire  collector  were  alu- 
minized and  overcoated  with  silicon  oxide  for  protection.  These  components  are  shown  in 
Fig.  30. 


Fig.  30.  Laser  head  components  for  hand-held  laser,  showing 
(from  top  to  bottom)  Xe  flash  lamp,  sapphire  pump  collector, 
and  NPP  rod. 


The  laser  head  was  mounted  in  a Vespel  holder  to  which  were  attached  the  l/2-in. -dia.  laser 
mirrors,  using  adjustable  O-ring  mounts.  The  present  laser  has  a semiconfocal  cavity  with  a 
13-cm-radius  high-reflecting  mirror  and  a flat,  10%  transmitting  output  mirror. 

The  laser-cavity  assembly  is  mounted  on  top  of  a power  supply  (2-1/2  x 3 x 5-3/4  in.)  con- 
taining a commercially  purchased  DC  inverter  with  LC  pulse-forming  network  and  SCR  trigger 
circuitry.  The  inverter  is  powered  by  rechargeable  12-V  Ni-Cd  batteries,  also  mounted  in  the 
power  supply,  capable  of  providing  one  thousand  1 -J  discharges  into  the  lamp.  The  laser  is 
fired  manually  by  a push-button  switch  at  rates  up  to  1 Hz.  The  complete  laser  is  shown  in 
Fig. 31. 


Fig.  31.  Complete  hand-held  NPP  laser 
with  cylindrical  laser  cavity  mounted  on 
top  of  power  supply. 


The  laser  rod  presently  in  use  is  1/2  x 1 x 7 mm.  The  internal  lamp  and  rod-mounting 
surfaces  on  the  sapphire  collector  have  a ground  finish,  which  we  have  not  index-matched  to 
the  lamp  and  rod.  Because  of  these  factors,  we  are  presently  limited  to  0.5  mj  output  with  1 J 
input.  Using  a 1-mm^  rod  and  better  cementing  techniques,  we  should  be  able  to  achieve  out- 
put energies  comparable  with  the  10-mJ  level  (at  1 J input)  obtained  with  our  laboratory -bench 
hollow-collector  laser  heads. 

The  size  of  the  present  resonator  assembly  is  determined  primarily  by  the  external  mir- 
rors and  mounts.  Using  miniature  mirrors  fastened  directly  to  the  sapphire  holder,  or  suitably 
polishing  and  coating  the  laser-rod  surfaces  themselves,  would  result  in  a much  smaller  laser 
assembly.  Significant  reduction  in  the  power-supply  size  should  also  be  achievable. 

3.  Q-Switched  Operation,  Stable  and  Unstable  Resonators 

The  next  phase  of  experimental  study  dealt  with  Q-switched  operation  of  the  NPP  lasers 
using  an  intracavity  saturable  absorber.  We  used  the  dye  bis-(4-diethylaminodithiobenzil)- 
nickel  (Kodak  14617)  having  a molecular  weight  of  685.68.  In  most  of  the  measurements,  the 
dye  was  dissolved  in  toluene  and  the  dye  solution  was  held  in  specially  fabricated  cells  inside 
the  laser  cavity.  This  was  done  primarily  for  the  convenience  of  being  able  to  vary  the  dye 
concentration  and  quickly  change  the  Q-switch  loss. 
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Fig.  32,  Spectrophotometer  absorption  curves  of  transmission  vs  wavelength  for 
2.5  x 10"5  M (solid  curve)  and  4.0  x 10-5  m (dashed  curve)  solutions  of  BDN  dye  in 
toluene.  The  path  length  is  1 mm  and  a 1-mm  cell  of  toluene  was  used  as  a reference. 
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A starting  1x10  M solution  of  dye  was  prepared  by  dissolving  14.37  mg  of  dye  powder 
(2.1  x 10  5 moles)  in  210  ml  of  toluene,  and  lower  concentration  solutions  were  then  made  by 
appropriate  dilution.  Figure  32  shows  two  spectrophotometer  absorption  curves  for  2.5  x 10"5M 
and  4.0  x 10  5 M solutions,  obtained  with  1-mm-thick  sample  cells  and  using  a reference  cell 
containing  pure  toluene. 


At  the  laser  wavelength  of  1052  nm.  the  respective  transmissions  are 

As  a check 

on  the  dilution  accuracy,  these  coefficients  have  a ratio  1.86/3.01  = 0.618,  compared  with 

0.625.  Using  the  absorption  coefficient 


0.83  and  0.74,  giving  absorption  coefficients  1.86  cm'1  and  3.01  cm  *,  respectively. 

i co 

the  concentration  ratio  2.5  x 10  5/4.0  x 10'5 

3.01  cm  1 for  the  4.0  x 10  5 M dye.  the  molecular  ground  state  absorption  cross  section  is 
calculated  to  be  aD  = 1.25  x 10  11>  cm2 


which  is  nearly  the  same  value  found  in  Ref.  14  for 

BDN  in  polymethylmethacrylate. 

Two  types  of  dye  cells  were  constructed.  The  first  was  made  using  two  l/2-in.-dia.  glass 
windows  separated  by  a thin  cylindrical  spacer  ring  1.4  mm  thick.  The  faces  of  the  spacer 
ring  were  polished  flat  and  parallel,  and  it  was  separated  along  a diameter  by  a slot  in  order 
to  accommodate  the  filling  and  venting  tubes.  The  windows,  spacer  sections,  and  tubes  were 
fastened  together  with  epoxy.  The  outer  faces  of  the  cell  windows  were  antireflection  coated, 
but  the  inner  faces  were  left  uncoated  because  the  refractive  index  mismatch  between  the  glass 
(~1.52)  and  the  toluene  solutions  (~1.50)  was  so  small.  This  cell  was  fastened  to  the  output 
mirror  holder  with  the  stainless -steel  fill  tube  at  the  bottom  of  the  cell  connected  by  Teflon 
tubing  to  a fill  syringe  containing  the  dye,  and  the  top  vent  tube  leading  to  a discharge  collection 
flask.  This  arrangement  allowed  initial  filling  of  the  cell  with  toluene  for  optimum  alignment, 
and  then  pumping  and  filling  with  dye  solutions  without  disturbing  the  optical  alignment. 

The  second  type  of  cell  was  identical  in  construction  to  the  first,  except  that  the  plane  inner 
surface  of  one  of  the  windows  was  dielectric  coated  to  form  the  output  mirror.  Two  cells  of 
this  type  were  made,  having  reflectivities  of  0.80  and  0.75  with  respect  to  the  toluene  filling 
medium. 
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Fig.  33.  Output  vs  input  energy  for  1-  x 
1-  x 7-mm  NPP  rod  in  10-cm  confocal 
cavity  with  20%  transmitting  output  mir- 
ror. The  open-circle  data  points  are  for 
the  normal -mode  cavity,  and  the  closed- 
circle  points  are  for  the  cavity  with  in- 
ternal type  I Q-switch  cell  containing 
pure  toluene. 


Fig.  34.  Relation  between  Q-switch  dye 
concentration  (left  scale)  and  threshold 
energy,  and  output  energy  (right  scale) 
vs  input  energy  for  single-pulse,  Q- 
s witched  output  inal-  xi-x  7-mm 
NPP  rod  with  confocal  cavity. 


The  initial  results  using  the  type  I dye  cell  with  the  1-  x 1-  x 7-mm  rod  in  the  confocal 
resonator  are  shown  in  Figs.  33  and  34.  The  effects  of  adding  the  dye  cell  itself  are  shown  in 
Fig.  33,  where  the  operation  with  and  without  the  toluene-filled  cell  are  shown  (the  normal- 
mode data  are  the  same  as  in  Fig.  22).  First,  the  data  of  Fig.  22  show  a threshold  increase  of 
about  100  mJ  per  10  percent  cavity  loss,  and  then  from  the  data  of  Fig.  33  we  estimate  an  added 
cell  loss  of  about  3.5  percent.  Note  that  this  is  too  small  to  be  evident  in  the  slope  efficiency, 
but  is  more  apparent  in  the  threshold  shift  with  introduction  of  the  dye  cell.  Next,  in  Fig.  34 
we  show  (open  circles,  left  scale)  the  threshold  as  a function  of  dye  concentration  in  the  cell. 
The  calculated  loss  coefficient  for  the  dye  is  0.75  cm-1  per  10*5  M,  or  for  the  present  cell, 

2 at  = 0.21  per  10  * M.  The  line  shown  through  the  data  has  a slope  of  200  mJ  per  10~5  M,  or 
about  0.20  loss  per  10  5 M,  although  the  fit  is  not  good,  and  differences  between  2 at  and  the 
round-trip  loss  1 -e  2ai  are  beginning  to  be  significant.  For  a given  concentration,  multiple 
Q -switched  pulses  with  separations  on  the  order  of  tens  of  microseconds  occurred  when  the 
input  energy  was  too  large.  In  measuring  the  single-pulse  Q-switched  energy  output,  we  op- 
erated with  an  energy  midway  between  the  threshold  and  double -pulsing  points  for  the  most 
stable  operation.  In  addition  to  the  energy  meter  at  the  laser  output  end,  we  put  a Si  diode  be- 
hind the  HR  mirror,  which  detected  a weak  output  leaked  through  this  mirror  and  was  used  to 
monitor  and  verify  single-pulse  operation  during  the  energy  measurements.  The  energy  of  the 
Q-swltched  pulses  measured  in  this  fashion  as  a function  of  pump  energy  is  also  plotted  in 
Fig.  34  (open  squares,  right  scale).  With  1 J input,  we  obtained  1.75  mJ  output  using  3 x lo'5  M 
dye.  This  represents  27  percent  of  the  6.4 -mJ  normal-mode  output  measured  with  the  cell 
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Fig.  3 5.  Time  dependence  of  Q-switched 
output  pulse  using  a 1-  x 1- x 7 -mm  NPP 
rod  and  confocal  cavity.  The  energy  is 
1.75  mJ,  and  pulse  width  is  ~7  nsec. 
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Fig.  36.  Relation  between  Q-switch  round-trip  loss  and  threshold  energy 
(bottom  scale)  and  single-pulse  output  energy  (top  scale)  for  a 1-  x 2-  x 
10 -mm  NPP  rod  in  a confocal  cavity. 


containing  toluene.  The  time  dependence  of  the  Q-switched  pulse  is  shown  in  Fig.  35,  as  mea- 
sured with  a fast  Si  p-i-n  photodiode  and  storage  oscilloscope.  The  width  of  the  pulse  (FWHM) 
is  ~7  nsec  (we  believe  the  small  second  pulse  is  a spurious  signal  probably  due  to  an  impedance 
mismatch). 

In  Fig.  36.  we  show  similar  data  obtained  with  the  1-  x 2-  x 10-mm  rod,  using  the  same 
dye  cell  and  confocal  resonator.  We  have  plotted  the  calculated  dye-cell  loss  coefficient  2 of 
(using  dye  concentrations  up  to  5 x 10  5 M)  on  the  vertical  axis,  and  the  solid  data  points  show 
the  energy  threshold  (bottom  scale)  for  varying  loss.  Even  though  we  are  in  a regime  where 
2 at  is  much  greater  than  1-e  2at  , and  even  exceeds  unity,  the  approximately  linear  relation 
is  found  because  the  gain  coefficient  2 g t which  balances  2 of  is  increasing  roughly  linearly  with 
input  energy.  In  other  words,  even  in  the  high-gain  regime,  the  round-trip  gain  e2at  when 

«2  fyl 

multiplied  by  e is  a constant  determined  by  the  other  fixed  cavity  losses.  As  before,  the 
Q-switched  output  energy  at  each  dye  concentration  was  measured  for  a pump  midway  between 
the  threshold  and  double-pulse  level,  but  keeping  EjN  less  than  1 J.  The  data  for  output  Q-switch 
energy  (open  circles,  top  axis)  are  shown  for  varying  dye  loss.  At  1 J input  with  5 x 10"5  M 
dye,  we  obtained  2.35  mJ  output. 

In  addition  to  the  liquid  cell,  we  also  used  some  plastic  sheet  Kodak  Q-switch  material 
(cellulose  acetate  containing  BDN  dye)  provided  by  J.  Gualtieri  of  USAERADCOM.  Although  no 
clear  plastic  sheet  was  available  as  a reference  for  the  spectrophotometer  transmission  ref- 
erence in  order  to  separate  intrinsic  absorption  loss  from  dielectric  reflection  loss,  we  esti- 
mate the  absorption  coefficients  to  be  in  the  range  given  by  the  data  points  with  brackets  in 
Fig.  36.  The  thresholds  and  output  energies  are  in  moderately  good  agreement  with  the  toluene- 
solvent  Q-switch  data,  although  we  believe  the  poorer  optical  quality  of  the  plastic  sheet  does 
result  in  larger  internal  laser  loss  and  lower  output.  Using  the  more  absorbing  sheet,  we  ob- 
tained 1.5  mJ  output  at  925  mJ  input. 

In  an  effort  to  obtain  single-mode,  low -divergence  output,  we  chose  to  use  the  type  of  cavity 
called  a hybrid  unstable  resonator.  Unstable  resonators  have  the  property  that  a geometric 
ray  trace  back  and  forth  through  the  laser  cavity  eventually  leads  to  the  ray  exiting  from  the 
cavity  around  the  edge  of  an  output  coupler.  The  most  common  types  of  such  resonators  have 
high-reflectance  mirrors  and  either  use  the  annular  output  from  around  one  of  the  end  mirrors 
or  insert  a scraper  mirror  within  the  cavity  to  deflect  this  annular  output.  The  ray-trace  con- 
cept is  equivalent  to  that  of  expanding  spherical  waves  which  have  the  desirable  feature  of  filling 
the  gain  medium  and  extracting  most  of  the  energy  within  the  cavity.  To  avoid  the  difficulties 
associated  with  fabricating  and  mounting  an  edge-coupled  output  mirror  or  introducing  an  addi- 
tional scraper  mirror,  the  hybrid  unstable  resonator  was  introduced  for  producing  uniform  low- 
divergence  output  in  small  Nd:YAG  laser  rangefinder  cartridges.17  In  this  design,  one  mirror 
is  partially  transmitting,  giving  a uniform  beam,  with  the  outer  portion  of  beam  corresponding 
to  the  conventional  annular  output  being  blocked  or  neglected. 

The  formulas  given  in  Ref.  17  for  estimating  the  magnification  and  geometric  loss  in  the 
geometric  (ray-tracing)  approximation  are 


and 
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where  we  have  assumed  a plane  output  mirror,  and  g = l t L/R,  with  L = cavity  length  and  R = 
convex  rear  mirror  radius  of  curvature.  With  the  10-  and  20-m-radius  mirrors,  the  calculated 
magnifications  and  edge  losses  are  M10  = 1.22  and  «10  = o.n.  M20  = 1.15  and  fi2Q  = 0.06. 

Data  for  the  hybrid  unstable  resonator  (HUR)  using  a 1-  x 2-  x 10-mm  rod  are  shown  in 
Fig.  37  for  normal-mode  (non-Q -switched)  operation.  Contrary  to  the  calculations,  the  20-m 
mirror  gave  a slightly  higher  threshold  and  less  energy  output  than  the  10-m  mirror,  even 
though  the  theoretical  geometric  loss  is  less  for  the  former.  This  may  be  related  to  differences 
in  the  equivalent-cavity  Fresnel  number,  which  is  closer  to  a half-integer  value  for  the  10-m 
mirror,  which  tends  to  reduce  the  geometric  loss.  The  experimental  results  with  the  10-m 
mirror  were  nearly  identical  with  either  the  21%  T output  mirror  or  the  20%  T type  II  Q -switch 
cell.  There  was  a significant  reduction  in  the  output  compared  with  that  of  the  stable  confocal 
resonator  (Fig.  23),  which  may  be  partly  due  to  the  geometric  loss  causing  higher  threshold  and 
lower  efficiency. 

The  data  for  Q -switched  operation  using  the  HUR  are  shown  in  Fig.  38,  presented  in  the 
same  manner  as  in  Fig.  36  for  the  stable  confocal  resorator.  With  an  input  energy  of  925  mj 
we  obtained  1.2  mj  Q-switched  output  using  5 x 10  5 M dye.  This  is  about  20  percent  of  the 
corresponding  normal-mode  HUR  output,  and  about  half  of  the  stable  resonator  Q-switched  output. 

Using  a larger  cross  section  2-  x 2-  x 10-mm  rod  and  greater  input  energies,  the  HUR  nor- 
mal mode  data  of  Fig.  39  were  obtained.  With  the  type  n Q cell,  18  mJ  output  occurred  at  2 J 
input,  compared  with  34.5  mJ  in  the  confocal  resonator  (Fig.  25).  When  Q-switched  using  the 
25%  T type  II  cell,  we  obtained  5.15  mJ  output  at  1.86  J input  (Fig.  40). 


Fig.  37.  Output  vs  input  energy  for  normal  - 
mode  operation  of  a 1 - x 2-  x 10-mm  NPP 
rod  in  a hybrid  unstable  resonator. 
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Fig.  38.  Relation  between  Q-switch  round- 
trip  loss  and  threshold  energy  (bottom  scale) 
and  single-pulse  output  energy  (top  scale)  for 
a 1-  x 2-  x 10-mm  NPP  rod  in  a hybrid  un- 
stable resonator. 


mode  operation  of  a 2-  x 2-  x 10-mm  NPP  trip  loss  and  threshold  energy  (bottom  scale) 
rod  in  a hybrid  unstable  resonator.  and  single-pulse  output  energy  (top  scale)  for 


a 2-  x 2-  x 10-mm  NPP  rod  in  a hybrid  un- 
stable resonator. 


Far-field  measurements  of  the  HUR  output  were  performed  using  two  techniques.  In  one. 
the  Scientech  power  meter  head  was  placed  6.4  m (21  ft)  from  the  laser  with  a variable  iris  in 
front,  and  "energy-in-the -bucket"  measurements  were  done,  with  the  divergence  angle  being 
calculated  from  the  aperture  opening.  The  other  method  used  an  infrared  display  measurement 
similar  to  that  for  the  near-field  fluorescence.  In  this  case,  a target  with  calibration  grid  was 
placed  at  the  far-field  iris  position  and  the  vidicon  camera  was  focused  on  the  target  and  inci- 
dent laser  spot.  The  image  of  the  spot  was  displayed  on  the  TV  monitor  and  scanned  with  the 
CVI  vtdeo  analyzer  to  give  a quantitative  profile  of  the  mode. 

In  principle,  one  should  use  a collimating  lens  at  the  laser  output  to  compensate  for  the 
geometric  divergence  of  the  beam,  which  has  a radius  of  curvature17 
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L + LR 


mirror  + z 


(42) 


where 


R = convex  mirror  radius 
z = distance  from  output  mirror 

At  z = 0,  the  calculated  radius  is  10  cm  x 1000  cm  = 1 m requiring  a 1 -m  lens  for  colllmation. 
We  did  not  have  this  long-focal-length  lens  available,  and,  in  any  event,  the  calculated  Fresnel 
diffraction  pattern  did  not  differ  significantly  with  or  without  the  lens  since  the  calculated  geo- 
metric divergence  is  small  compared  with  diffraction. 
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RELATIVE  INTENSITY 


A typical  far-field  scan  is  shown  in  Fig.  41.  This  particular  pattern  was  obtained  with 
normal-mode  operation  of  the  1-  x 2-  x 10-mm  rod  in  the  HUR  using  the  10-m  mirror  and 
type  II  toluene -filled  dye  cell  with  20%  transmission.  The  spatial  calibration  was  found  from  a 
separate  CVI  scan  of  the  illuminated  grid  pattern  on  the  screen.  Most  of  the  energy  falls  within 
a 2.4-mrad  angle,  or  a 1.5-cm  spot  diameter.  The  calculated  spot  size  for  a 1-mm  source  is 
1.2  cm  for  a diffraction -limited  beam.  The  observed  spot  does  not  exhibit  the  expected  diffrac- 
tion lobes,  although  these  should  be  less  than  5 percent  of  the  peak  amplitude,  and  also  appears 
to  be  symmetric,  even  though  the  rod  cross  section  is  rectangular.  Using  the  power  meter  for 
this  same  measurement  gave  energies  of  6.0  mj  in  the  near  field  and  4.8  mJ  and  5.3  mJ  in  2.3- 
and  2.7-mrad  apertures,  respectively.  Under  the  same  conditions  which  previously  gave  1.05  mJ 
Q-switched  output  in  the  near  field,  we  obtained  0.75  mJ  in  3 mrad  in  the  far  field,  with  a scan 
pattern  almost  identical  to  Fig.  41.  This  Q-switched  output  was  more  erratic  in  that  the  spot 
shape  or  mode  pattern  was  occasionally  quite  asymmetric,  although  its  general  size  remained 
the  same. 

The  same  measurements  were  repeated  for  the  2-  X 2-  x 10-mm  AR  rod  (near-field  data  of 
Fig.  40).  At  the  beginning  of  the  measurements,  alignment  for  maximum  near-field  energy  gave 
5.5  mj  Q-switched  output  at  1.92  J in.  However,  the  far  field  showed  a double-lobed  pattern. 
When  the  alignment  was  adjusted  to  give  a symmetric  far-field  output,  we  obtained  3.9  mj  in  a 
3.1 -mrad  aperture.  The  CVI  far-field  scan  of  this  output  is  shown  at  the  top  of  Fig.  42,  and  the 
time  behavior  is  shown  at  the  bottom  (6  nsec  FWHM). 


IV.  SUMMARY  AND  CONCLUSIONS 


As  the  first  part  of  this  study,  we  have  experimentally  measured  the  anisotropic  thermal 
conductivity  and  specific  heat  of  pure  NdP5C>14.  With  these  parameters,  in  principle,  it  is 
possible  to  calculate  the  time-  and  space -dependent  temperature  distribution  of  NPP  laser  rods, 
given  the  heating  function  found  from  the  pump  absorption  and  the  boundary  conditions  for  cool- 
ing the  rod.  We  may  not  be  able  to  specify  these  quantities  precisely,  because  of  effects  such 
as  convective  heating  and  cooling  of  the  rod  and  the  complex  geometry  of  the  pump  absorption, 
but  nevertheless,  simple  approximations  should  give  good  semi-quantitative  predictions  of  the 
heating  to  be  expected.  The  next  important  step  should  be  to  calculate  and  measure  how  such 
heating  affects  both  the  lasing  performance  and  the  crystal  integrity  of  the  laser  rod. 

The  lasing  performance  is  ffected  by  temperature  from  both  intrinsic  fluorescence  prop- 
erties and  macroscopic  optical  properties.  The  two  dominant  processes  affecting  the  intrinsic 
ionic  fluorescence  .'roper ties  are  the  increase  in  thermal  population  of  the  *1^/2  man‘f°W-  re" 
quiring  a larger  population  inversion,  and  the  increased  linewidth  of  the  lasing  transition 

at  higher  temperatures,  lowering  the  effective  laser  cross  section  and  gain.  For  a uniform 
temperature  rise,  the  first  factor  would  cause  a practical  change  in  the  resonant  loss  term  of 
(AE/kT)  (l/T)  which  at  room  temperature  is  approximately  3%/°C.  Since  the  round-trip  reso- 
nant loss  itself  is  about  5%/cm,  a temperature  rise  on  the  order  of  20°C  would  be  necer^ary 
to  make  a significant  change  in  laser  performance,  since  the  output  coupling  should  far  exceed 
the  resonant  loss.  The  transitions  in  NPP  are  predominantly  homogeneously  broadened  from 
a phonon  Raman  scattering  process.  The  theoretical  linewidth  is  proportional  to  T2,  giving  a 
fractional  decrease  in  the  laser  cross  section  or  gain  of  2/T  or  0.7%/*C  at  room  temperature. 
This  slight  decrease  should  be  offset  to  some  extent  from  the  increased  overlap  of  adjacent 
transitions  which  contribute  to  the  gain  at  a given  wavelength  near  the  R^-Y^  dominant  peak. 

If  operation  of  the  laser  is  required  at  temperature  exceeding  50’C,  it  may  be  necessary  to 
incorporate  temperature -compensation  circuitry  or  other  means  to  avoid  or  allow  for  multiple 
Q-switched  pulsing  which  could  occur  at  the  low -temperature  extreme,  where  the  laser  gain 
would  be  higher  and  loss  lower. 

The  other  direct  thermal  effects  on  laser  performance  arise  from  changes  in  the  macro- 
scopic optical  properties  such  as  the  index  of  refraction.  In  Nd:YAG,  which  because  of  its 
cubic  symmetry  possesses  no  intrinsic  birefringence,  thermal  gradients  lead  to  thermally  in- 
duced birefringence  which  limits  the  mode  quality  and  output  power.  In  NPP,  the  intrinsic 
birefringence  is  significant,  so  that  the  polarization  properties  are  unlikely  to  be  affected  by 
temperature  changes.  However,  the  large  thermal  gradients  may  cause  changes  in  the  highest 
gain  index  component,  n^,  large  enough  to  lead  to  thermal  focusing  or  defocusing  effects  in  the 
rod,  changing  the  effective  resonator  parameters  and  output  divergence.  Some  initial  measure- 
ments of  the  change  in  optical  path  length  with  temperature  have  been  reported,  but  the  individual 
contributions  of  dn/dT  and  df/dT  have  not  been  determined. 

This  latter  property  is  also  important  in  determining  the  mechanical  behavior  of  the  laser 
rod  as  its  temperature  changes.  To  avoid  straining  the  rod,  the  thermal  expansion  coefficients 
of  the  mount  should  match  those  of  the  rod.  In  view  of  the  large  anisotropy  of  the  thermal  con- 
ductivity. some  degree  of  anisotropy  in  the  thermal  expansion  should  be  expected  and  the  mount 
designed  accordingly  if  it  is  to  operate  over  a wide  temperature  range.  In  addition  to  the  pos- 
sible strains  produced  at  the  rod -mount  Interface,  it  may  also  be  necessary  to  examine  the 


stresses  induced  within  the  rod  from  a combination  of  thermal  gradients  and  anisotropic  thermal 
expansion  under  extreme  operating  conditions. 

Our  ab  initio  calculations  of  the  time-dependent  NPP  laser  output  give  good  qualitative 
agreement  with  the  experimental  results.  The  quantitative  differences  arise  from  the  unknown 
parameters  of  the  system,  primarily  the  exact  absorption  geometry  and  the  transverse  laser 
mode  distribution.  Because  the  absorption  coefficients  are  large  in  NPP,  both  of  these  factors 
are  critical  in  determining  the  laser  behavior,  as  demonstrated  by  the  calculated  results  for  a 
slightly  reduced  mode  size.  Improvements  to  the  calculations  could  be  made  by  more  accurate 
modeling  of  the  time-dependent  lamp  radiation  and  by  calculating  the  rate-equation  solutions  in 
the  high-gain  regime,  but  in  view  of  the  previously  mentioned  unknown  quantities,  this  would 
probably  serve  little  purpose. 

Experiments  with  the  flash-lamp-excited  NPP  lasers  have  demonstrated  the  capability  of 
such  devices  to  provide  several  millijoules  of  Q-switched  energy  in  a beam  with  less  than 
4 mrad  divergence.  These  parameters  do  not  fully  demonstrate  the  capability  of  such  lasers, 
since  many  aspects  of  the  laser  operation  were  not  optimized.  In  order  to  summarize  the  re- 
sults with  a view  toward  extrapolating  the  final  device  potential,  we  shall  review  some  of  the 
separate  aspects  of  the  laser  performance.  First,  a minimum  value  of  normal-mode  energy 
conversion  of  2.1  percent  total  (2.5  percent  differential)  has  been  demonstrated  with  a 
2-  x 2-  x 20-mm  rod  operating  at  3. 3 J input,  70  mJ  output.  This  result  was  obtained  in  a 
confocal  resonator  aligned  for  uniform,  multimode  output  with -an  uncoated  rod.  Other  mea- 
surements have  shown  no  differences  between  uncoated  rods  and  rods  that  had  anti-reflection 
coatings.  These  were  equivalent  in  the  sense  that  coating  did  not  reduce  any  intracavity  losses, 
as  determined  from  the  laser  threshold  and  output  efficiency.  However,  from  the  behavior  of  the 
2-  x 2-  x 20-mm  rod  at  the  highest  pump  levels,  it  appears  that  anti -reflection  coating  may  be 
necessary  to  suppress  laser  oscillation  from  an  uncoated  facet.  With  such  coating,  an  output 
of  88  mJ  at  4 J input  should  be  attainable. 

This  output  might  be  improved  by  redesign  of  the  laser  head  to  provide  more  uniform  illu- 
mination of  the  laser  pod  and  a more  uniform  gain  distribution  within  the  rod.  Techniques  such 
as  use  of  elliptical  focusing  might  prove  effective,  at  the  cost  of  a slightly  larger  head  cross 
section.  Although  our  initial  experiments  with  dual-lamp  excitation  were  inconclusive,  proper 
design  of  such  a laser  head  should  also  be  effective  in  yielding  more  uniform  inversion. 

Another  area  of  improvement  is  in  the  design  of  the  hybrid  unstable  resonator.  A compar- 
ison of  the  data  in  Figs.  26  and  39.  for  example,  shows  a differential  efficiency  of  2 percent  for 
the  confocal  resonator  and  1.5  percent  for  the  hybrid  unstable  resonator,  with  the  HUH  having 
60  percent  of  the  confocal  output.  This  H1JR  output  is  a near-field  measurement,  and,  typically, 
85  percent  of  this  is  collected  in  a 3-mrad  aperture  in  the  far  field,  reducing  the  slope  efficiency 
and  relative  output  to  1.3  and  50  percent,  respectively.  If  one  uses  the  calculated  geometric 
loss  factor  of  11  percent,  the  theoretical  maximum  relative  efficiency  should  be  T/L  + T = 0.20/ 
(0.11  + 0.20)  = 65  percent.  Different  resonator  parameters  or  an  improved  design,  possibly 
using  a conventional  unstable  resonator  with  a small-reflecting-area  output  coupler,  should  be 
investigated.  It  should  be  feasible  to  extract  -75  percent  of  the  multimode  output  in  a diffraction- 
limited  beam,  for  an  estimated  66  mJ  output  at  4 J input  in  normal-mode  operation. 

The  major  remaining  problem  is  extraction  of  the  normal-mode,  long-pulse  energy  in  a 
short  Q-switched  pulse.  This  problem  has  been  examined  in  a more  general  analytic  fashion 
by  considering  the  losses  from  the  Q-switch  in  the  initial  unbleached  state  and  above  threshold 
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in  the  saturated  state.  Our  approach  has  been  somewhat  more  empirical,  but  more  quantitative 
since  these  two  dye  loss  terms  have  been  treated  as  dynamic  variables  in  the  laser  rate  equations, 
and  the  integrated  effect  of  the  time-varying  loss  has  been  found.  The  particular  data  of  Fig.  21 
show  that  roughly  0.4  to  0.5  of  the  normal-mode  output  should  be  achievable  for  the  present  type 
of  saturable -absorber  dye.  In  practice  (Figs.  39  and  40),  this  ratio  is  about  5 mj/18  mJ  = 0.28 
(near  field).  One  factor  which  could  cause  a reduction  in  Q-switch  efficiency  is  the  nonuniformity 
of  the  laser  mode,  which  could  prevent  full  dye  saturation  at  the  mode  periphery.  Assuming 
this  Q-switch  efficiency  factor  to  be  0.3,  a predicted  energy  output  is  then  0.3  x 66  mJ  = 20  mJ 
at  4 J input.  Therefore,  for  this  goal  of  20  mJ  in  a diffraction-limited  Q-switch  pulse  to  be 
reached,  one  must  have  a minimum  reduction  factor  of  ~0.25  of  the  experimentally  measured 
80-mJ  long-pulse  multimode  energy.  In  our  experiments,  this  overall  measured  reduction 
factor  is  ~4  mJ/35  mJ  = 0.11,  so  an  improvement  of  about  40  percent  is  necessary  in  both  the 
HUR  and  Q-switch  efficiencies. 

Effects  of  the  dye  cells  on  the  mode  divergence  have  not  been  measured.  Nonlinear  index 
changes  related  to  the  dye  bleaching  may  act  as  an  internal  lens,  affecting  the  far-field  output, 
and  should  be  considered  in  choosing  the  dye  cell  parameters  of  length  and  concentration. 

Another  influence  on  these  parameters  is  the  possibility  of  optical  damage.  In  one  case  where 
a semi -concentric  cavity  was  tested,  the  small  mode  area  at  the  dye  cell  position  caused  visible 
damage  to  the  cell.  Examination  under  a microscope  showed  that  the  damage  occurred  at  the 
inner,  uncoated  window-dye  interface  and  presumably  resulted  from  destructive  absorption  by 
the  dye,  with  possible  combined  effects  of  high  thermal  loading  and  chemical  decomposition 
causing  damage  spots  on  the  glass  surface.  While  this  was  an  extreme  case,  with  a power 
density  much  higher  than  would  be  reached  in  a non-focused  mode,  it  indicates  the  possibility 
of  infrared  damage  to  the  dye  if  energy  absorption  density  is  too  large.  Such  a situation  might 
be  reached  for  example  in  thin  cells  with  high  dye  concentration,  even  for  non-focused  modes. 

Another  necessary  step  is  the  improvement  of  the  optical  quality  of  the  plastic  sheet  dye 
Q-switches.  Experimentally,  they  have  yielded  only  0.50  to  0.7  5 of  the  energies  obtained  with 
comparable  density  liquid  Q-switch  cells.  Likely  causes  of  this  reduction  are  poorer  surface 
flatness  and  parallelism  and  material  inhomogeneities  in  the  plastic  sheet.  Since  the  plastic 
sheet  we  used  was  not  AR  coated,  the  non-parallelism  and  misalignment  of  the  sheet  may  also 
have  introduced  internal  losses  to  the  laser. 

Although  our  calculations  have  shown  that  the  maximum  Nd  concentration  is  required  for 
the  most  efficient  lasing  performance,  this  result  may  not  be  true  with  combined  effects  of 
reduced  mode  size  and  intensity -dependent  Q-switch  loss.  The  slight  loss  in  overall  absorption 
efficiency  with  lower  Nd  concentration  may  be  offset  by  the  increased  gain  uniformity  in  the 
laser  rod  and  its  effect  on  the  mode  profile.  Unfortunately,  the  calculation  of  an  unstable  res- 
onator mode  with  non-uniform  gain  is  quite  complex,  and  at  the  present,  there  are  no  (Nd,  La)PP 
laser  rods  for  obtaining  experimental  data. 
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APPENDIX 

GREEN'S  FUNCTION  SOLUTION  OF  THE  ANISOTROPIC  HEAT  FLOW  EQUATION 


The  solution  of  the  heat-flow  equation  (1)  can  be  found  by  integrating  the  solution  to  the 
Green's  function  equation  (2)  over  the  time-  and  space -dependent  heating  source  function.  The 
solution  to  (2)  is  found  as  follows.  Let  the  spatial  Fourier  transform  of  G equal  T(p*).  Then 


(AxPx  + Vy  + Vz  * F " pC  IT  = 4*6(T) 

where  r = t - t . The  solution  of  Eq.  (Al)  is 


_ . . 4»  I 'AA  + Vy  + Azpz2! 

r = u(t)  ^ exp ft* r 


(Al) 


<A2) 


The  solution  G is  the  inverse  transform  of  T, 
G(r. 


r.r).  If  C eip'rr(r.T)dr 

(2*)  J-« 


(A3) 


where,  for  example,  the  contribution  to  the  factorable  integral  from  the  x component  is 


“ ipxX  *AvPv  T/pC 

e x e x x 


-pCx 


Using  similar  results  for  the  y and  z coordinates. 


j0pr/xr{- &(*:*£,*  k)\  ■ 


<A4) 


(A5) 


This  solution  is  the  impulse  response  to  a point-source  excitation  at  t = t — t = 0 so  that  for  a 
sinusoidal  driving  function  we  must  perform  the  convolution  over  cos  wto> 


p1  Pn 

T<r-t)sJ  ?7  CObwtoG(t~to) 
With  the  change  of  variable 


dt_ 


(A6) 


the  integration  over  { gives 
T(r 


P , 2 2 2 .-1/2  f / 

•t)3  47  (t  + f + A")  (AxAyV*  / R*  exp  f iwt  - ./ISpC 


(A7) 


This  is  the  solution  for  a time-periodic  point-source  distribution,  which  must  be  integrated 
to  find  the  response  from  the  present  line  source.  In  our  geometry,  the  line  source  lies  in  the 


plane  at  z = 0,  and  extends  in  the  y direction.  The  temperature  sensor  is  located  a distance 
"d"  from  the  heater  source,  and  is  also  in  the  plane  at  z = 0.  Since  "d"  is  very  much  less  than 
the  heater  length,  we  will  neglect  edge  effects,  and  take  the  limits  of  y integration  from  -«  to 
+«.  The  line-source  response  is 


Using  the  change  of  variable 


riu>pC 


from  the  integral  definition  of  the  modified  Bessel  function,  Kq. 

The  actual  heater  and  sensors  are  not  lines  of  zero  thickness,  but  have  some  finite  width  w 
When  this  factor  is  integrated  as  in  Ref.  5,  an  additional  factor  of  [sinh(aw/2)/(aw/2)|^  is  found 
using  arguments  as  above,  where  the  tactor  n now  has  the  anisotropic  tensor  element  of  the 
thermal  conductivity. 


a = (iuipC/A  » 


and  the  voltage  response  of  the  sensor  is 


(A12) 


with  all  quantities  defined  as  in  Eq.  (6).  Although  both  the  magnitude  and  phase  of  V( 
sured,  it  suffices  to  use  only  the  magnitude  in  determining  the  unknowns  a and  A.: 


W ( I ® I w/2r  J ° 

The  numerical  calculation  of  KQ  is  found  using  the  series  expansions  of  Ref.  18 


K0(«)  ■ -|l°g(*/2)  + yj  IQ(*)  + £ 

m=l 


where  the  modified  Bessel  function  lQ(z)  can  be  found  from 

2/>  ,m 


° (mir 


m=0 


and 


y = Euler's  constant  = 0.57721...  . (A15) 

In  our  calculations,  the  complex  argument  z is  of  the  form  n/T  v,  with  v real,  so  we  will  de- 
termine the  real  and  imaginary  parts  of  the  Bessel  functions  with  this  choice  of  variable.  From 
Eq.  (A14)  we  define 


S(z) 


2 ate?  d.i ....  .a. 

(ml)2  2 m 


(A16 ) 


m=l 


giving 


. * 2 / < >2  m tl  . « >m  » > 

S(z)  = i l lziL_  II*  * ' I 

(2m  + l!)2  2 2m  + 1 


m=0 


2 ivW'(-.r  ,±, 


( A17) 


= i Im  S + Re  S 


and 


V (v2/4)2m  (-l)m  , . v2  y (v2/4)2m  (~l)m 

_ I?  4 “ 


IQ(z)  = 


m=0 


(2m!  )“ 


1=0  (2m  + 1!) 


(A18) 


= Re  I + i Im  I 
o o 

For  computational  purposes  we  define  the  quantities 
1 


f = 


(ml) 


vZ/4  = y 


cm  = ~LZ  (T  + 7 + • • ' + m> 

m (mi)2  * 2 m 

2 

y = x 


so  that  the  above  functions  can  be  put  into  the  following  forms  for  efficient  numerical  analysis: 
IQ(*)  = 1 - x(f2  - x(f4  - x(ffc. ..)))  + iy  [fj  - x(f3  - x(f5  - x(. . . )))] 

S(z)  = -x(c2  - x(c4 ))  + iy  [Cj -x  (c3 -x(c& -x(...  )))J  . (A19) 

where  sufficient  terms  are  taken  to  insure  convergence.  Then, 

Kq(z)  = -(log  (v/2)  + y)  Re  IQ  + | Im  I + Re  S 


(A20) 
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